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FOREWORD 


The  theory  for  predicting  the  motions  of  a  ship  in  a  seaway  is  one  of  the 
triumphs  of  research  in  ship  hydrodynamics.  Given  a  rather  small  amount  of 
information  about  a  ship  and  the  seaway,  one  car.  predict  heave  and  pitch  motions 
to  a  remarkable  degree  of  accuracy  without  recourse  to  model  tests  or  empirical 
data.  Lateral-plane  motions ,  sway  and  yaw,  can  also  be  predicted  with  reasonable 
accuracy. 

However,  when  one  tries  to  predict  roll  motion,  one  realizes  what  good  luck 
we  have  had  in  analyzing  heave,  pitch,  sway,  and  yaw:  These  are  not  sensitive 
to  the  effects  of  fluid  viscosity.  Roll  motion  is  extremely  sensitive  to  viscosity 
effects,  especially  to  viscosity  -  induced  flow  separations.  In  addition,  roll 
motion  is  strongly  influenced  by  the  presence  of  bilge  keels,  which  are  difficult 
to  analyze  even  by  the  classical  methods  of  hydrodynamics  of  an  ideal  fluid. 

During  the  1978-79  academic  year,  our  Department  was  fortunate  in  having  as 
a  visiting  scholar  Professor  V.  Himeno  of  the  University  of  Osaka  Prefecture.  He 
is  well-known  in  Japan  for  his  research  on  viscous-fluid  problems  of  ship  hydro¬ 
dynamics,  and  the  Department  of  which  he  is  a  member  is  distinguished  for  its 
research  on  ship  roll  damping. 

Therefore  I  was  especially  pleased  when  he  agreed  to  my  request  that  he 
prepare  a  report  describing  the  state  of  the  art  in  predicting  roll  damping. 

As  in  many  areas  of  naval  architecture,  Japan  is  in  the  forefront  in  developing 
practical  procedures  for  predicting  ship  roll  damping.  Professor  Himeno  has  been 
closely  associated  with  these  developments. 

As  he  makes  clear  in  this  report  there  are  many  aspects  of  this  problem 
that  have  not  yet  been  adequately  analyzed.  However,  in  the  great,  tradition  of 
Japanese  naval  architecture  research,  theory  is  used  as  far  as  possible,  and  the 
gaps  are  filled  with  empirical  information.  More  research  is  needed,  but  a 
usable  procedure  for  predicting  roll  damping  is  describe:  . 

In  the  appendix,  a  computer  program  is  presented  for  predicting  roll  damping. 
This  program  from  Osaka  Prefecture  University  was  tested  at  The  University  of 
Michigan  by  having  an  undergraduate  compile  it  and  use  it.  The  information  pro¬ 
vided  in  the  Appendix,  together  with  the  comments  built  into  the  program,  were 
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sufficient  for  this  student  to  use  the  program. 


T.  Francis  Ogilvie 
Department  of  Naval  Architecture 
and  Marine  Engineering 
The  University  of  Michigan 


-vi- 


ACKNOWLEDGMENT 


NOMENCLATURE 


[Note:  Numbers  in  parentheses  indicate  equations  where  more  information 
can  be  found  about  the  quantity  listed.] 

Ar  Wave -amplitude  ratio  (Coe  Fig.  4.6) 

A<j,  Inertia  coefficient  in  roll  equation  of  motion  (2.1) 
a.  Extinction  coefficient  (2.15) 

B  Ship  beam 

Be  Damping  coefficient  component:  bare-hull  eddy  making  (4.1) 

Bp  Damping  coefficient  component:  bare-hull  skin  friction  (4.1) 

Bl  Damping  coefficient  component:  lift  (due  to  forward  speed)  (4.1) 

Bv*  Damping  coefficient  component:  bare-hull  wave  making  (4.1) 

Bek  Damping  coefficient  tompone.  „ :  total  of  bilge-keel  pressure 
effects  (4.2) 

Beeh  Damping  coefficient  component:  hull  pressure  change  due  to  bilge 
keels  (4.1) 

Bern  Damping  coefficient  component:  bilge-keel  normal  pressure  (4.1) 
BBkw  Damping  coefficient  component:  bilge-keel  wave  making  (4.1) 
Equivalent  linear  damping  coefficient  (2.5) 

B<j,  Roll-damping  term  in  equation  of  motion  (2.1) 

Bj  Coefficients  in  expansion  of  ,  j  =  1,2,...  (2.2) 

b  Extinction  coefficient  (2.15) 

bBK  Width  of  bilge  keel 

bi  Effective  width  of  bilge  keel  (4.20) 

Cg  Hull  block  coefficient 

Cq  Drag  coefficient  (4.17)  (See  Fig.  4.9) 

Cp  Pressure-difference  coefficient  (4.7) 
cj  Cp  in  front  of  bilge  keel 

Cp  Cp  behind  bilge  keel 

C,j,  Restoring- force  coefficient  in  equation  of  motion  (2.1) 

c  Extinction  coefficient  (2.15) 

d  Ship  draft 

Fn  Froude  number 

f  Empirical  coefficient  giving  velocity  increment  at  the 
bilge  circle  (4.20) 

<34  Metacentric  height  (restoring-moment  lever  arru)  (2.1) 
g  Acceleration  of  gravity 


KG 

Distance  from  keel  to  center  of  gravity 

Of 

ship 

k 

Reduced  frequency,  eL/U 

L 

Ship  length 

^BK 

Length  of  bilge  keel 

Roll  excitation  moment  (2.1) 

N 

"N-coefficient , "  (2.19) 

N10 

Value  of  N  for  $A  =  10°  (2.20) 

0 

Origin  of  coordinates 

OG 

Distance  downward  from  origin  to  center 

of 

gravity 

(4.5) 

R 

Bilge  radius  of  hull 

r 

Mean  distance  from  center  of  gravity  to 

bilge  keels 

(4.5) 

rmax 

Maximum  distance  from  roll  axis  to  hull 

surface  (4. 

7) 

So 

Width  of  distribution  of  Cp  on  hull 

Tn 

Natural  period  of  roll  (2.4) 

t 

Time  variable 

U 

Forward  speed  of  ship 

w 

Weight  of  ship 

a 

B1/2A(j)  (2.4) 

ae 

Be/2A^  (2.9) 

a 

B2/h*  (2.4) 

Y 

B3/A4  (2.4) 

Aip 

*n-l-+n  (2'1Sa) 

Radiation  wave  amplitude  (Fig.  4.6) 

Scale  ratio  of  ship  to  model 

\j 

Kinematic  viscosity  of  water 

Sd 

(u?d/g  (4.15) 

P 

Density  of  Water 

a 

Area  coefficient  of  a  cross-section  of  the 

hull 

T 

Oai/g  (4.15) 

4> 

Roll  angle  (2.1) 

<f>A 

Amplitude  of  roll  motion  (2.1) 

to 

Frequency  (rad/sec) 

Special  Notations  : 

V  Displacement  volume  of  ship 

Indicates  nondimens ion al  form  of  quantity, 
[subscript]  indicates  value  at  zero  forward  speed 
indicates  2-D  value  for  a  cross-section  of  hull. 
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1. 


INTRODUCTION 


Roll  motion  is  one  of  the  most  important  responses  of  a  ship  in  waves. 

The  roll  motion  of  a  ship  can  be  determined  by  analysing  various  kinds  of 
moments  acting  on  the  ship,  virtual  and  actual  mass  moments  of  inertia,  roll 
damping  moment,  restoring  moment,  wave  excitation  and  other  momemts  caused  by 
other  modes  of  ship  motion.  Among  them,  the  roll  damping  moment  hari  been  con¬ 
sidered  to  be  the  most  important  term  that  should  be  correctly  predicted.  It 
is  needed  not  only  at  the  initial  stage  of  ship  design  to  secure  the  safety 
of  a  ship,  but  also  to  obtain  a  better  understanding  of  ship  motions  in  waves. 

Since  the  age  of  W.  Froude ,  a  number  of  theoretical  and  experimental 
works  has  been  made  concerning  the  predictions  of  roll  damping  and  roll  mo- 
tion  of  ships.  The  recent  development  of  the  "Strip  Method"  has  made  it  pos¬ 
sible  to  calculate  almost  all  the  terms  in  the  equations  of  ship  motions  in 
waves  with  practical  accuracy,  except  for  t.he  roll  damping.  The  necessity 
of  obtaining  the  roll  damping  of  ships  has  been  pointed  out  in  the  recent 
recommendations  of  the  Seakeeping  Committee  of  the  International  Towing  Tank 
Conference  (ITTC) -  Notwithstanding  these  efforts,  it  seems  that  a  complete 
solution  of  this  problem  has  not  yet  been  reached. 

Difficulties  in  predicting  the  roll  damping  of  ships  arise  from  its 
nonlinear  characteristics  (due  t.o  the  effect  of  fluid  viscosity)  as  well  as 
from  its  strong  dependence  on  the  forward  speed  of  ship.  Moreover,  the  fact 
that  these  various  effects  have  influences  on  the  value  of  roll  damping  that 
are  of  the  same  order  of  magnitude  makes  the  problem  even  more  complicated 
in  the  absence  of  bilge  keels. 

After  the  classical,  works  by  Bryan  [  1  ]  and  Gawn  (  2  ]  ,  we  can  recog¬ 
nize  an  epochmaking  period  a  couple  of  decades  ago  in  the  history  of  research 
on  roll  damping.  Experiments  on  bilge  keels  by  Martin  f  3  ]  ,  Tanaka  [  4  ] , 
and  Kato  [  5  ) ,  theoretical  works  on  the  vortex  flow  near  bilge  keels  by 
S&sajiina  (  6  ] ,  consideration  of  hull-friction  damping  by  Kato  [  7  ] ,  and 
study  of  the  surface-tension  effect  by  Ueno  [  8  ] ,  all  of  these  works  ap¬ 
pear  in  this  period. 

Furthermore,  we  can  cite  here  Hishida's  theoretical  studies  [  9  )  on  the 
wavemakir.g  roll  damping  due  to  hull  and  bilge  keels  and  Hanaoka's  mathematical 
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formulation  [10]  for  the  wave  system  created  by  an  oscillatory  motion  of 
an  immersed  flat-plate  wing  with  low  aspect  ratio.  Also,  an  extensive 
series  of  free-roll  tests  at  zero  ship  speed  for  ordinary  ship  hull  forms 
was  carried  out  by  Watanabe  and  Inoue  [11]  and  tests  at  forward  speed  by 
Yamanouehi  [12] „  Some  of  the  results  of  these  works  have  often  been  used 
even  up  to  the  present  time  or  have  given  background  to  recent  works.  This 
fact  cannot  but  remind  us  again  of  the  difficulty  of  treating  the  ship  roll¬ 
damping  problem  rigorously. 

It  can  be  said  that  the  recent  works  started  about  a  decade  ago,  mainly 
associated  with  the  experimental  check  of  the  accuracy  of  the  strip  method. 
Much  data  on  radiation  forces  acting  on  ship  hull,  including  roll  damping, 
have  been  accumulated  through  the  forced  oscillation  tests  carried  out  by 
Vugts  [13] ,  Fujii  and  Takahashi  [14] ,  Takaki  and  Tasai  [15] ,  and  Takezawa 
et  a 2.  [16] . 

These  experiments  have  clarified  that  there  are  still  considerable  dif¬ 
ferences  between  measured  values  of  roll  damping  and  those  predicted  by  ex¬ 
isting  methods.  In  this  period  much  effort  has  also  been  made  f  GIT  Cbtaliiiii^ 
snip  roll  damping,  for  example,  works  by  Bolton  [17] ,  Lofft  [18] ,  and 
Lugovski  et  al.  [19]  ,  concerning  the  effect  of  bilge  keels,  Gersten's  studies 
[20]  on  the  viscous  effects,  and  free-roll  experiments  hy  Takaishi  et  al. 
t21]  and  Tanaka  et  al.  [22]. 

Moreover,  what  should  be  noted  here  is  the  extensive  and  systematic  works 
in  Japan  that  have  recently  been  carried  out  through  the  cooperation  of  the 
Japan  Shipbuilding  Research  Association,  especially  in  the  Committees  of 
SR108,  succeeded  by  SR125,  131  and  161  [23].  In  the  prediction  method  of 
ship  roll  damping  considered  there,  damping  is  divided  into  several  com¬ 
ponents,  for  instance,  friction,  eddy,  lift,  wave,  and  bilge-keel  con^ponents. 
Then  the  total  damping  is  obtained  by  summing  up  these  component  dampings  pre¬ 
dicted  separately.  This  attempt  appears  to  have  had  a  certain  success  for 
ordinary  ship  hull  forms. 

The  objective  of  this  article  is  to  describe  the  present  state  of  the 
art  in  these  recent  attempts  as  well  as  other  existing  formulas  for  ship  roll 
damping.  Furthermore,  for  convenience  in  ship  design,  it  is  intended  that 
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the  available  expressions  and  formulas  should  be  described  in  full  detail  as 
much  as  possible,  so  that  their  values  can  be  calculated  promptly  once  the 
particulars  of  a  ship  are  given. 

In  Chapter  2  the  various  methods  of  representation  for  roll  damping 
coefficients  and  their  relationships  are  stated  and  rearranged  in  terms  of  an 
equivalent  linear  damping  coefficient.  Then,  for  prediction  methods  of  roll 
damping,  simple  methods  are  introduced  in  Chapter  3,  including  the  use  of  data 
from  a  regression  analysis  of  model  experiments.  In  Chapter  4,  the  newest 
treatment  for  component  dampings  is  stated,  and  available  formulas  for  each 
component  are  fully  described  there.  Comparisons  are  made  of  measured  and  pre¬ 
dicted  total  damping,  which  is  the  sum  of  the  component  dampings. 

Chapter  5  concerns  the  prediction  methods  for  ship  roll  motion.  However, 
it  is  not  the  full  present  state  of  the  art.  The  description  is  limited  to 
the  problem  of  how  to  use  the  formulas  of  nonlinear  roll  damping  in  order  to 
obtain  the  solution  of  the  roll  equation  of  motion  in  regular  or  irregular 
sea.  Finally  an  example  of  a  FOKTBAN  statement  of  a  computer  program  for  the 
component  dampings  is  given  in  the  Appendix. 


2. 


REPRESENTATION  OP  ROLL  DAMPING  COEFFICIENTS 


Many  ways  of  representing  roll  damping  coefficients  have  been  used,  de¬ 
pending  on  whether  roll  damping  is  expressed  as  a  linear  or  nonlinear  form, 
which  form  of  the  non-dimensional  expressions  is  to  be  used,  and  by  what 
experimental  method  its  value  was  measured,  for  instance,  forced-roll  test 
or  free-roll  experiment.  Some  of  the  expressions  most  coi'monly  used  are 
introduced  here  and  the  relationships  among  them  are  reviewed  and  rearranged 
in  terms  of  a  linearized  damping  coefficient. 

2.1  Nonlinear  Damping  Coef ficients 

The  equation  of  roll  motion  has  recently  been  expressed  as  a  three-degree- 
of-freedom  form, including  sway  and  yaw  motions  simultaneously.  However,  in 
order  to  limit  the  discussion  here  to  the  problem  of  nonlinear  roll  damping, 
we  can  write  down  the  equation  of  the  roll  motion  of  a  ship  in  the  following 
simple  single-degree-of-freedom  form: 

A^i  +  (i )  +  =  M^Wt)  ,  (2.1) 

In  Eq. (2.1),  g  represents  the  toll  angle  (with  the  amplitude  <$^) ,  A^  the 
virtual  mass  moment  of  inertia  along  a  longitudinal  axis  through  the  center 
of  gravity  and  the  coefficient  of  restoring  moment,  which  is  generally 

equal  to  W*GM  (W  the  displacejhent  weight  of  the  ship  and  GM  the  metacentric 
height) .  Furthermore,  tfy  stands  for  the  exciting  moment  due  to  waves  or 
external  forces  acting  on  the  ship,  w  the  radian  frequency  and  t  the  time. 
Finally,  B^  denotes  the  roll  damping  moment,  which  is  now  considered. 


Although  only  the  main  terms  of  roll  motion  have  been  taken  into  account 
in  Eq.  (2.1,,  coupling  terms  being  neglected,  it  can  be  said  that  Eq.  (2.1) 
almost  corresponds  to  that  of  three  degrees  of  freedom  when  we  consider  the 
wave  excitation  term  as  the  Froude  form  with  a  coefficient  of  effective 

wave  slope.  This  is  because  the  concept  of  the  effective  decrease  of  wave 
slope  turns  out,  after  Tasai's  analysis  [24],  to  correspond  to  the  effect  of 
the  sway  coupling  terms. 


We  can  express  the  damping  moment  as  a  series  expansion  of 

and  1 4>  |  in  the  form 


-5- 


3^  =  B14>  +  B2iU!  +B3$2h-...  ,  (2.2) 

which  is  a  nonlinear  representation.  The  coefficients  B^  ,  B2  ,  . . .  in 
Eq.  (2.2)  are  considered  as  constants  during  the  motion  concerned.  In  other 
words,  these  values  may  possibly  depend  on  the  scale  and  the  mode  of  the 
motion,  for  instance,  on  the  amplitude  $A  and  the  frequency  a)  when  the 
ship  is  in  a  steady  roll  oscillation. 

Dividing  Eq. (2.1)  by  ,  we  can  obtain  another  expression  per  unit 

mass  moment  of  inertia: 

<f>  +  2cuf)  +  $ | $  +  -,43  +  =  m.  (<ot)  ,  (2.3) 


where 


n 


m<j> 


(2.4) 


In  Eq.  (2.4) ,  the  quantities  «n  and  Tn  represent  the  natural  frequency 
and  the  period  of  roll,  respectively. 


A  term  of  the  form  |/|||  might  be  added  to  the  right-hand  side  of 
Eq. (2.2).  This  term  corresponds  to  the  effect  of  surface  tension  at  water 
level  of  the  ship  hull.  Ueno  (  8  ]  investigated  this  effect  and  concluded 
that  the  surface  tension  might  cause  a  considerable  error  in  the  values  of 
the  damping  coefficients  when  a  small  model  is  used  in  small  amplitude,  of 
oscillation.  However,  this  effect  is  not  considered  hereafter,  because  the 
surface  tension  depends  strongly  on  the  condition  of  the  painted  surface  of 
the  model  hull  as  well  as  on  that  of  the  water  surface,  and  because  it  can 
be  neglected  in  the  case  of  roll  amplitude  with  moderate  magnitude  for  a 
ship  model  of  ordinary  size. 

To  obtain  the  values  of  these  coefficients  of  nonlinear  damping  directly 
through  a  steady-state  forced-roll  experiment,  in  which  <f>A  and  u>  are 
specified,  we  would  probably  need  numerical  techniques  to  fit  the  solution 
of  the  assumed  equation  to  the  measured  data.  Such  an  attempt  does  not  seem 
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to  have  been  done.  Instead,  the  usual  way  that  has  been  taken  is  to  assume 
some  additional  relations  concerning,  say,  energy  consumption ,  linearity  of 
damping  and  its  independence  of  <f>  ,  which  will  be  described  later. 

2.2  Equivalent  Linear  Damping  Coefficients 

Since  it  is  difficult  to  analyse  strictly  the  nonlinear  equation  stated 
in  the  preceding  section,  the  nonlinear  damping  is  usually  replaced  by  a 
certain  kind  of  linearized  damping  as  follows : 

B^(4>)  =  Be*  .  (2.5) 

The  coefficient  Be  denotes  the  equivalent  linear  damping  coefficient. 
Although  the  value  of  depends  in  general  on  the  amplitude  and  the 

frequency,  because  the  damping  is  usually  nonlinear,  we  assume  Be  is 
constant  during  the  specific  motion  concerned. 

There  are  several  ways  to  express  the  coefficient  B  in  terms  of  the 

e 

nonlinear  damping  coefficients  B^_  ,  B2  and  so  on.  The  most  general  way 

is  to  assume  that  the  energy  loss  due  to  damping  during  a  half  cycle  of  roll 
is  the  same  when  nonlinear  and  linear  dampings  are  used  [24] .  If  the  motion 
is  simple  harmonic  at  radian  frequency  u j  , 

Be  =  Bl+^-w,fAB2  +|-w24>|b3  .  (2.6) 

For  more  general  periodic  motion,  Eq.  (2.6)  can  be  derived  by  equating  the 
first  terms  of  the  Fourier  expansions  of  Eqs.  (2.5)  and  (2.2)  [15]. 


lor  convenience  in  analyzing  the  equations  of  lateral  motions,  the 
nondimensional  forms  of  these  coefficients  are  defined  as  follows: 


for  i=lr2,3 


(2.7) 


where  p  ,  V  and  B  stand  for  fluid  density,  displacement  volume  and 

breadth  of  ship,  respectively.  Then  Eq.  (2.6)  takes  the  following  nondimen¬ 
sional  form: 
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8  " 

+  3^T^AB2 


+ 


^AB3 


(2.8) 


Corresponding  to  Eq. (2.3),  we  can  define  an  equivalent  linear  damping 
coefficient  ae  =  Ee/2A ^  per  unit  mass  moment  of  inertia: 


=  ct  +  ~-m<K6  + 


(2.9) 


Since  these  coefficients  still  have  dimensional  values  (except  for  the  second, 
6),  the  following  dimensionless  forms  are  often  used,  especially  for  the 
linear  terms  ae  and  a  . 


(2.10) 


In  case,  of  irregular  roll  motion,  there  is  another  approach  to  lineariza¬ 
tion  of  the  roll  damping  expression.  After  the  works  of  Kaplan  (26 ]  and 
others  [27]  [23] ,  we  assume  that  the  difference  of  the  damping  moment  between 

its  linearized  and  nonlinear  forms  can  be  minimized  in  the  sense  of  the  least 
squares  method.  Neglecting  the  term  B3  for  simplicity,  we  define  the  discrepancy 
<5  in  the  form  : 


6  =  +B2|$|£ -Be$  .  (2.11) 

Thru  we  can  minimize  E{62},  the  ertpectation  value  of  the  square  of  6 
during  the  irregular  roll  motion,  assuming  that  the  undulation  of  the  roll 
angular  velocity,  4>  ,  is  subject  to  a  Gaussian  process  and  that  the  coeffi¬ 

cients  Be  ,  B1  and  B2  remain  constant: 

3|-  E{62}  =  -2(B1-Be)E{|2)  -  2B2E{j>2  |$|  }  =  0  .  (2.12) 

After  some  calculations  we  can  reach  the  form 

B„  »  B,  +-CT.B  ,  (2.13) 

e  1  *  4  2 

where  the  factor  o.  represents  the  variance  of  the  angular  velocity  ip 
It  is  claimed  in  recent  works  [23]  that  this  form  is  useful  for  analyzing 
roll  behavior  in  irregular  seas. 


Moreover,  as  an  unusual  way  of  linearization,  we  can  equate  the  nonlinear 
expression  to  the  linear  one  at  the  instant  when  the  roll  angular  velocity 
takes  its  maximum  value  during  steady  oscillation: 

Be  =  Bj  +0)4  aB  .  (2.14) 

This  form  seems  to  correspond  to  a  collocation  method  in  a  curve-fitting 
problem,  whereas  Eq„ (2.6)  corresponds  to  the  Galerkin  approach.  Since  there 
is  a  difference  of  about  15%  between  the  second  terms  of  the  right  hand  sides 
of  Eqs .  (2.6)  and  (2.14),  the  latter  form  may  not  be  ■'’■alia  for  the  analysis 
of  roll  motion.  But  it  may  be  used  as  a  simple  way  of  analyzing  forced-oscil¬ 
lation  test  data  to  obtain  the  values  of  these  coefficients  promptly  from 
the  time  history  of  the  measured  roll  moment. 

However,  the  most  common  way  to  obtain  these  nonlinear  damping  coeffi¬ 
cients  through  forced  oscillation  tests  is,  first,  to  find  the  equivalent 
linear  coefficient  Be  in  Eq. (2.5)  by  assuming  that  the  forced-roll-test 
system  is  subject  to  a  linear  equation,  and,  second,  to  fit  Eq. (2.6)  to  the 
Be  values  obtained  by  several  test  sequences  with  the  amplitude  <1>A  varied. 
Then  we  can  obtain  the  values  of  these  damping  coefficients,  ,  B2  and 

so  on,  which  are  independent  of  the  amplitude  of  roll  oscillation. 

It  should  be  noted  here  that  this  condition,  the  independence  of 
amplitude,  was  not  stated  when  we  derived  Eq. (2.6).  Therefore  we  might 
obtain  different  values  of  the  coefficients  from  the  original  ones  if  the 
coefficients,  especially  ,  should  depend  on  the  amplitude,  particularly 

in  the  presence  of  bilge  keels.  We  should  keep  these  things  in  mind  when  we 
use  a  formula  like  Eq. (2.6) . 

2.3  Extinction  Coefficients 

A  free-roll  test  is  probably  the  simplest  way  to  measure  roll  damping 
of  ship  or  model.  In  a  model  test,  sway  and  yaw  motions  are  usually  restrained 
to  avoid  the  effect  of  the  horizontal  motions.  On  the  other  hand,  heave 
and  pitch  motions  are  often  kept  free  to  avoid  the  error  due  to  the  sinkage 
force  in  the  presence  of  forward  speed,  although  it  is  of  course  desirable  to 
make  the  vertical  motions  as  small  as  possible.  The  roll  axis  is  usually 
taken  through  the  center  of  gravity  of  the  model,  the  radius  of  gyration  of 
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the  model  is  adjusted  to  the  value  of  the  actual  ship  considered,  and  the 
restoring  moment  lever  GM  is  also  measured  through  a  static  inclination 
test. 

In  a  free-roll  test,  the  model  is  rolled  to  a  chosen  angle  and  then 
released.  The  subsequent  motion  is  measured.  Denote  by  the  absolute 

value  of  roll  angle  at  the  time  of  the  n-th  extreme  value.  The  so-called 
extinction  curve  expresses  the  decrease  of  <$n  as  a  function  of  mean  roll 
angle.  Following  Froude  and  Baker,  we  fit  the  extinction  curve  by  a  third- 
degree  polynomial : 

A<p  =  a<j>m  +  b<f)^  +  ctf)^  (deg.)  ,  (2.15) 


where 

Ad>  =  4n_^  -  <j>n  ,  (2.15a) 

=  Wn-1  +*n]/2  *  (2.15b) 

The  angles  are  usually  measured  in  degrees  in  this  process. 

The  coefficients  ,  b  and  c  are  called  extinction  coefficients.  The 
relation  between  these  coefficients  and  the  damping  coefficients  can  be 
derived  by  integrating  Eq. (2.1)  without  the  external- force  term  over  the  time 
period  for  a  half  roll  cycle  and  then  equating  the  energy  loss  due  to  damping 
to  the  work  done  by  restoring  moment.  The  result  can  be  expressed  in  the 
form 

_  (l)_  /  q  -5  1 

A<j>  =  I B  + -%  $_B  +4i)2lj)2B  (rad)  .  (2.16) 

Y  2  i  377  n’a  2  4  nYm  3J 

Comparing  Eq.  (2.16)  with  Eq.(2.15)  term  by  term,  we  can  obtain  the  relations 
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It  should  be  noted  here  again  that  the  condition  for  the  validity  of  Eq. 
(2.17)  is  that  the  coefficients  ,  3?  ,  and  a  ,  8  ,  ....  should  be  inde¬ 

pendent  of  the  roll  amplitude.  As  we  can  see  in  the  later  chapters,  the 
effect  of  bilge  keels  appears  mainly  in  the  term  B,  ,  and,  further,  the 
value  of  varies  with  roll  amplitude.  In  such  a  case,  Eq. (2.17)  will  not 

remain  valid.  Only  the  part  of  B,  which  is  independent  of  the  amplitude 
is  related  to  the  coefficient  b-  The  other  part  of  B2  that  is  inversely 
proportional  to  the  anplitude  will  apparently  be  transferred  tc  the  coeffi¬ 
cient  a  ,  and  the  part  proportional  to  the  amplitude  will  appear  in  c.  In 
place  of  a  term-by-term  comparison,  therefore,  it  will  probably  be  reasonable 
to  define  an  equivalent  extinction  coefficient  ae  and  to  compare  it  with 
the  equivalent  linear  damping  coefficient  Bg  as  in  the  form 

*e  “  a  +  b*m  +  c*m  =  >r  Be  •  <2-18) 

$ 

We  are  also  familiar  with  Bertin' s  expression  [25],  which  can  be  writ¬ 
ten  in  the  form 


-  N<{!2  (deg.)  .  (2.19) 

The  coefficient  N  can  be  taken  as  a  kind  of  equivalent  nonlinear  expression, 
and  it  has  been  called  an  "N-coefficient. "  As  seen  from  Eq. (2.15) , 

N  =  ^ — +  b  +  cd  (deg.)  .  (2.20) 

"b  m 

The  value  of  N  depends  strongly  on  the  mean  roll  angle  <j>m  ,  so  that  its 
expression  is  always  associated  with  the  value,  being  denoted  as  , 

N2q  and  so  on,  where  value  of  N  when  <J>m  *  10°,  etc. 


3.  PREDICTION  OF  ROLL  DAMPING:  I  .  SIMPLE  METHOD 


When  the  principal  dimensions  of  a  ship  form  are  given,  the  most  reli¬ 
able  way  to  obtain  the  roll  damping  of  the  ship  at  present  time  seems  to  be 
to  carry  oat  a  model  experiment.  Since  the  scale  effect  of  the  damping  is 
considered  to  be  associated  mainly  with  the  skin  friction  cn  the  hull,  which 
makes  a  small  contribution  to  total  damping,  the  data  from  the  model  tests 
can  easily  be  transferred  co  the  actual  ship  case  by  using  an  appropriate 
nondimensional  form  of  roll  damping,  for  instance,  Eq. (2.8). 

If  model-test  data  are  not  available,  it  is  necessary  to  estimate  the 
roll-damping  value  by  using  certain  kinds  of  prediction  formula.  There  are 
two  different  ways  of  estimation  at  present  time.  One  is  to  obtain  an  empiri¬ 
cal,  experimental  formula  directly  through  the  analysis  of  model  tests  on 
actual  ship  forms.  The  other  is  to  break  down  roll  damping  into  several 
components  and  then  estimate  the  value  by  summing  up  the  values  of  those 
components  individually  predicted. 

The  latter  is  considered  to  be  more  rational,  so  that  it  has  become  the 
recent  trend  of  approach  in  Japan.  To  begin  with,  however,  some  examples  of 
the  former  approach  are  described  in  this  chapter  in  order  to  know  the  mag¬ 
nitude  of  ship  roll  damping  easily. 

3.1  Watanabe-Inoue-Takahashi  Formula 

A  couple  of  decades  ago,  Watanabe  and  Inoue  [11]  established  a  formula 
for  predicting  the  roll  damping  of  ordinary  ship-hull  forms  at  zero  advance 
speed  in  normal-load  condition,  on  the  basis  of  both  an  extensive  series  of 
model  tests  and  some  theoretical  considerations  on  the  pressure  distribution 
on  the  hull  caused  by  ship  roll  motion.  Their  original  formula  has  been 
modified  slightly  by  them  [28]  so  as  to  be  applicable  to  a  wider  range  of 
ship  forms,  including  ships  with  large  values  of  block  coefficient. 

Takahashi  [29]  proposed  a  form  of  forward-speed  modification  multiplier 
to  be  applied  to  the  value  at  zero  ship  speed,  thus  expressing  the  advance- 
speed  effect  on  roll  damping.  We  may  call  this  approach  the  Watanabe-Inoue- 
Takahashi  formula. 


-1-1- 


-12- 


This  can  be  expressed  in  terms  of  an  equivalent  linear  damping  coeffi¬ 
cient  of  the  form 

wn 

Be  =  Be0  [1  +  0 . 8{  1-exp (l-10Fn)  }  —  ]  f  (3.1) 

(a2 

where  Beo  stands  for  the  value  of  Be  at  zero  ship  speed.  Its  value  can 
be  expressed  in  terms  of  the  extinction  coefficients  a  and  b  ,  as  follows: 


®eo  ~  2“eA<t 


2  ui 

~  “nA<f>U  +: 


“n"<fi v “  ($^:deg)  .  (3.2) 

Furthermore  these  coefficients  can  be  related  to  the  values  of  the  N- 
coefficier.ts  N, 


(3.3) 


10 

and 

N20 

at  the 

roll  . 

a 

+  t), 

a 

10 

=  10 

2 

to 

o 

II 

20+b 

and  the  virtual  mass  moment  of  inertia  A^  can  be  calculated  by  the  strip 
method  or  otherwise  determined  through  the  approximate  relationship  of  the 
natural  roll  frequency  to  both  the  coefficient  A^  and  the  restoring  moment 
of  the  ship. 

The  N-coeff icients  can  be  expressed  after  Watanabe  and  Inoue  in  the 

form 


'10 


20 


+ 


n 

ao 


1.5 


1.0 


(3.4) 


where  L,B  and  d  represent  ship  length,  beam  and  draft.  The  quantity 
Ajj  denotes  the  area  of  the  bilge  keel  at  one  side  of  the  hull.  The  distance 
l  is  defined: 

%  =  KG  -  |  ,  (3.5) 

and  the  quantity  f  is  a  function  of  waterline  area  coefficient  c  : 

w 

f  =  1_  JL  +_s _ 4__  +^_. 

m+1  2m+l  3m+l  4mtl  ' 


'  w 


1-Cw 


(3.6) 
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The  first  square  bracket  an  the  right-hand  side  of  Eq , (3.4)  corresponds 
to  a  kind  of  drag  coefficient  of  nonlinear  damping.  The  coefficient 
and  n,, 0  stand  for  those  of  the  naked  hull  at  tile  roll  amplitudes  10° 
and  20°  ,  expressed  in  the  form 


)  1  0.78  ) 

CBd 

^  0.03  ) 

i  ■  L ! 

L 

+  )  C . 02  i 

f 

(3.7) 

where  Cg  represents  ship  block  coefficient.  Finally,  the  quantity 
corresponds  to  a  certain  efficiency  cf  the  bilge  keels,  the  value  of  which 
can  be  determined  from  Fig.  3.1,  given  by  the  authors  themselves  as  a  function 
of  c  and  the  aspect  ratio  b  JL  of  the  biloe  keel. 

B  BIv/  Bis. 

We  car.  thus  obtain  the  magnitude  of  the  roll  damping  easily  once  the 
principal  dimensions  of  a  ship  and  its  bilge  keels  are  given-  Since  these 
formulas  were  established  on  the  basis  of  a  rather  large  amount  of  experimental, 
data,  it  can  be  hoped  that  they  will  offer  a  reasonable  estimation  of  roll 
damping  in  trie  early  stage  of  ship  design.  As  an  example,  the  comparison  of 
T?k abash i  's  formula  with  some  experimental  data  is  shown  ir.  Fig.  3.2.  The 
agreement  seems  to  be  acceptable  in  the  ordinary  range  of  ship  speed. 

o 


Figure  3.2.  Effect  of  advance  speed  on  roll  damping  force 
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Tiie  background  of  these  formulas  should  also  be  mentioned.  Eq.  (3.2)  can 
be  derived  from  Eq. (2.17),  which  represents  the  term-by-term  comparison  between 
the  coefficients  Bj  ,  B2  »  ...  and  the  extinction  coefficients  a  ,  b 
and  c  .  For  its  validity,  however,  the  as sumption sate  needed  that  the 
coefficients  Bj  and  B2  take  their  values  at  end  are  constant,  that 

is,  are  independent  of  frequency  and  are  not  affected  by  the  roll  amplitude 
-j>A  .  These  assumptions,  especially  the  former,  might  cause  a  misprediction  for 
the  ship  rolling  in  a  frequency  range  away  from  the  ship  natural  frequency. 

Therefore  these  formulas  should  be  applied  to  the  case  of  a  normally 
loaded  ship,  and  then  only  near  the  natural  frequency,  where  ship  rolling 
usually  becomes  important. 

3.2  Tasai-Takaki ' s  Table 

As  a  second  example  of  a  simple  method,  we  cite  the  results  of  Tasai  and 
Takaki 1  s  experiment  [23],  which  has  recently  been  carried  out  for  the  purpose 
of  obtaining  typical  values  of  roll  damping  for  ordinary  ship  hull  forms.  The 
roll  damping  of  four  typical  kinds  of  ship  form  was  measured  by  forced-roll 
test  at  specified  values  of  both  Froude  number  and  the  roll  amplitude,  with 
frequency  varied.  The  data  obtained  were  fitted  by  regression  analysis  to 
the  forms  of  Eq. (2.8),  including  two-  and  three-term  damping  coefficients. 


Table  3.1  Particulars  of  Models 


container 

cargo  ship 

ore  carrier 

tanker 

length  Lpp  (m) 

3.00 

3.00 

3.00 

3.00 

breadth  B  (m) 

0.4354 

0.4783 

0.493 

0.4719 

draft  d  (m) 

0.1628 

0.1957 

0.194 

0.1828 

disp^jic^ment 

121.61 

199.84 

233.4 

220.57 

0.5717 

0.7119 

0.8243 

0.8519 

cM 

0.970 

0.9905 

0.9975 

0.9946 

GM  (m) 

0.017 

0.02174 

0.050 

0.06077 

OB  (m) 

-0.0425 

-0.025 

0.089 

0.0993 

Kj/Lpp 

0.239 

0.2172 

0.2356 

0.2494 

Kb/B 

0.382 

0.3240 

0.2602 

0;2513 

^BK/^pp 

0.25 

0.25 

0.25 

0.37 

b3K/B 

0.0148 

0.0159 

0.0142 

0.00869 
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Table  3.1  shows  the  principal  dimensions  of  the  models,  and  Tables  3.2 
and  3.3  present  the  nondimens ional  values  of  the  nonlinear  damping  coeffi¬ 
cients,  taken  up  to  the  second  and  third  terms,  respectively 


Table  3.2  Damping  coefficients  of  2nd-order  approximation 


As  an  example  for  the  validity  of  these  curve  fittings,  the  comparison 
with  the  experiment  on  the  C&  =  0.71  cargo  ship  model  is  shown  in  Fig.  3.3. 

The  agreement  seems  to  be  quite  reasonable,  even  in  the  region  of  comparatively 
large  roll  amplitude.  Since  the  result  concerns  a  typical  ship  form,  it 
should  be  quite  useful  frr  the  prediction  of  magnitude  of  ship  roll  damping 
at  the  initial  stage  of  ship  design.  It  can  be  applied  to  other  cases  of 
different  ship  form  and  ship  speed  by  interpolating  or  extrapolating  the  values 
in  the  tables. 

What  should  also  be  noted  here  is  the  limitation  on  the  application  of 
these  tables.  It  is  assumed  that  the  damping  coefficients  B1  ,  and 

in  Eq. (2.8)  are  constant  for  specified  ship  form  anu  Froude  number.  There¬ 
fore  the  tables  do  not  cover  the  case  without  bilge  keels,  in  which  wave 
damping  might  prevail  and  the  magnitude  of  the  nonlinear  terms  might  not  be 
so  large. 
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Table  3.3  Damping  coefficients  of  3rd-order  approximation 


ISH 

ore  carrier 

tanker 

container 

cargo  ship 

B, 

0  .00  308 

0  .00209 

0  .00082 

0  .0006  1 

0 

E2 

0  .03262 

0  .04168 

0  .0  3690 

0  .04908 

b3 

0  .12170 

0  .0  3877 

0  .084  74 

0  .08994 

0  .00  359 

0  .00316 

0  . 10 

'  Kj|  «■ 

0  .04110 

0  .04453 

0  .07783 

0  .0  3581 

A 

Bi 

0  -00  344 

0  .00  374 

0 .00242 

0.15 

b2 

0 .04254 

0.02531 

0  -03755 

b3 

0 .05524 

0 .09835 

0 .08755 

Bi 

0  .00332 

0.20 

b2 

0  .0  3551 

l 

b3 

0.05226 

( 

ih 

0.00628 

0.00  309 

!  0.25 

j  ~ 

B2 

0.02125 

0.04033 

b3 

0.03567 

0.02206 

0.00671 

0.275 

B2 

0.01402 

B3 

0.05097 

We  have  seen  a  couple  of  examples  of  simple  methods  for  predicting  the 
order  of  magnitude  of  ship  roll  damping.  Other  similar  methods  might  also  be 
available  for  that  purpose.  To  improve  the  accuracy  of  the  prediction  method, 
however,  we  should  consider  the  phenomenon  of  roll  damping  from  a  much  more 
physical  and  hydrodynamic  point  of  view.  In  this  sense,  the  concept  of 
component  damping,  which  will  be  described  in  the  next  chapter,  will  present 
more  fruitful  results  applicable  to  a  wider  variety  of  cases. 


3rd  order  approximation 

2nd  order  approximation 


Figure  3,3. 


4.  PREDICTION  OF  ROLL  DAMPING: 


XI.  COMPONENT  ANALYSIS 


4.1  Definition  of  Component  Damping 

It  has  been  pointed  out  since  a  couple  of  decades  ago,  as  mentioned  pre¬ 
viously,  that  roll  damping  of  a  ship  is  caused  by  various  kinds  of  fluid  flow 
phenomena,  for  instance,  skin  friction  of  the  hull,  eddy  shedding  from  the 
hull,  free-surface  waves,  etc.  It  has  also  been  noted  that  roll  damping  is 
strongly  affected  by  the  presence  of  bilge  -keels,  rudder  and  appendages. 

In  recent  days,  the  concepts  of  these  effects  have  been  considerably 
clarified  and  research  work  on  the  individual  effects  has  been  carried  out. 

The  objective  of  this  chapter  is  to  define  all  these  effects  from  the  recent 
point  of  view  of  the  concept  of  component  dampings  and  to  describe  available 
prediction  formulas  for  them. 

To  begin  with,  we  assume  that  the  total  roll-damping  coefficient  for  an 
ordinary  ship  hull  form  can  be  divided  into  seven  components,  that  is,  friction, 
eddy,  lift  and  wave  damping  for  naked  hull,  and  normal-force  damping  of  bilge 
keels,  hull-pressure  damping  due  to  bilge  keels,  and  wave  damping  of  bilge 
keels,  Thi3  assun^tion  can  be  expressed  in  terms  of  equivalent  linear  damping 
coefficients  of  the  form 

Be  =  BF  *  Be  +  SL  +  B.W  +  Bbkn  +  3bkh  +  BBkw  ,  (4.1) 


or  we  can  sum  up  the  bilge  keel  terms, 


Be  =  %  +  %  J'  bl  +  %  +  bBK  •  (4'2) 

Although  these  coefficients  are  seemingly  linear,  their  values  may  vary 
with  the  roll  amplitude  and  the  frequency  to  .  For  brevity  the  effect 

of  appendages,  except  for  rudder  and  bilge  keels,  is  not  considered  here. 

The  rudder  is  assumed  to  be  included  in  the  main  hull  configuration. 

We  can  define  these  component  dampings  as  follows,  neglecting  or  including 
their  mutual  Interactions: 

Friction  damping  B_  is  caused  by  the  skin-friction  stress  on  the  hull 
£ 
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in  roll  motion,  so  that  it  may  possibly  be  influenced  by  the  presence  of  waves 
and  bilge  keels. 

The  eddy  damping  BE  stands  for  the  nonlinear  damping  (  B2 1  $  |  <p  ) 
caused  by  the  pressure  variation  on  the  naked  hull,  excluding  the  effect  of 
waves  and  bilge  keels.  However  it  may  apparently  include  the  linear  term 
B  ,  which  would  be  transferred  from  the  nonlinear  term  if  a  part  of  B2 

were  inversely  proportional  to  the  amplitude  $  .  In  the  presence  of  ship 

forward  speed,  eddy  damping  represents  the  nonlinear  part  of  the  lifting  effect 
of  the  hull  itself  in  roll  motion,  whereas  the  linear  part  is  defined  as  the 
lift  damping  B  . 

ij 

'The  wave  damping  B^  denotes  the  increment  of  the  hull-pressure  damping 
due  to  the  presence  of  free  surface  waves,  so  that  it  includes  the  interaction 
between  waves  and  eddies  and  between  waves  and  lift.  However,  since  these 
interactions,  will  be  very  small,  it  will  be  almost  linear. 

Bilge  keel  damping  B  represents  the  increment  of  pressure  damping  due 
tc  the  presence  of  a  pair  of  bilge  keels.  This  term  consists  of  the  following 
three  components:  The  normal- farce  damping  of  bilge  keels  B0KN  is  due  to 
the  normal  force  on  the  bilge  keels  themselves.  The  second  is  hull-pressure 
damping  due  to  bilge  keels  B  ,  which  corresponds  to  the  pressure  change 
on  the  hull  when  bilge  keels  are  installed.  Therefore  this  term  stands  for 
an  interaction  between  hull  and  bilge  keels.  The  rest  is  the  wave  damping  of 
bilge  keels  3DVTI  .  Since  the  first  two  terms  do  not  concern  the  free  surface, 
this  term  represents  the  change  of  values  of  bbkn  BBKH  <iue  to  waves- 

This  term  also  includes  interaction  between  the  hull  (with  bilge  keels)  and 
waves- 

We  have  presented  so  many  component  dampings  that  it  is  natural  to  ask 
which  component  is  the  most  important.  The  answer  is  quite  difficult  because , 
as  will  be  seen  later,  almost  all  the  components  (except  Bp  and  )  41:0 

of  the  same  order  of  magnitude  for  an  ordinary  ship  huli  form.  This  is  the 
very  fact  that  has  made  the  prediction  of  ship  roll  damping  difficult. 

The  subdivision  of  roll  damping,  as  stated  above,  may  not  always  be  based 
on  the  hydrodynamic  point  of  view  but  may  stand  rather  on  a  practical  basis 
for  convenience  in  predicting  the  roll  damping  and  carrying  out  the  ■'uariments. 
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However,  we  can  still  note  here  some  remarks  on  these  components.  The  terms 
Bj.  ,  Bw  and  B  w  can  safely  be  treated  as  nonviscous  damping  ,  while  the 
others  can  be  regarded  as  viscous  damping. 

Tangential  Stress  on  the  hull  contributes  only  to  the  term  Bp  .  The 
rest  of  the  damping  .is  caused  by  the  normal  stress,  that  is  ,  the  pressure  on 
the  hull  or  on  the  bilge  keels.  The  surface-wave  effects  appear  mainly  in  the 
terms  B,,  and  ■  The  other  terms  can  be  considered  as  free  from  waves , 

although  the  term  Bp  also  includes  a  small  wave  effect  in  its  definition. 

To  distinguish  linear  damping  from  nonlinear  damping  is  quite  difficult 
at  this  stage.  Of  course,  the  non-viscous  dampings  like  and  Bw  can  be 

regarded  as  linear,  but  also  some  parts  in  the  terms  B„  and  might 

be  linear,  due  to  their  dependence  on  Reynolds  number  or  the  Keulegan-Carpenter 
number.  The  forward-speed  effect  is  included  in  all  terms. 

These  characteristics  will  be  much  better  clarified  after  the  individual 
prediction  methods  are  described  in  the  following  sections. 

Hereafter  the  subscript  o  represents  the  case  of  zero  speed,  the  super- 
script  the  value  in  the  cross-sectional  plane.  The  roll  axis  passes  through 
the  center  of  gravity  unless  otherwise  stated. 

4.2  Friction  Damping 

In  predicting  the  value  of  friction  damping,  we  ignore  the  effect  of 
waves  and  regard  the  ship  hull  form  as  an  equivalent  axisymmetric  body,  for 
which  the  dimensions  will  be  defined  later.  Then  the  skin  friction  laws  for 
a  flat  plate  in  steady  flow  are  applied  to  the  roll  motion  of  the  body,  as 
shown  by  Kate  [  7  ] ,  Takaki  [30] ,  and  Schmitke  [31] . 


We  cite  Kata's  formula  as  an  example.  In  the  absence  of  forward  speed, 
Kato  applied  Blasius'  formula  for  laminar  flow  and  Hughes'  formula  for  turbu¬ 
lent  flow  to  the  peripheral  boundary  layer  of  a  circular  cylinder.  Kato’s 
formula  can  be  expressed  in  terms  of  an  equivalent  linear  damping  coefficient 
as  follows : 


Bycf  =  0.  787pSr|/w\>  |l  +  0.00814  — 


•  306 


(<]>A  :rad. ) 


(4.3) 
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where  p  and  v  represent  density  and  kinematic  viscosity  of  fluid.  In  Eq. 
(4.3),  the  first  term  in  the  brackets  gives  the  result  for  the  case  of  laminar 
flow,  which  is  used  for  the  naked  model  hull,  while  the  second  term  gives 
the  modification  for  turbulent  flow  by  Hughes'  formula,  applicable  to  both  the 
model  hull  with  bilge  keels  and  the  actual  ship  hull.  The  quantities  S  and 
rg  were  originally  defined  as  the  surface  area  and  the  radius  of  the  cylinder. 
For  the  case  of  a  ship  hull,  however,  they  represent  the  wetted  surface  area 
and  the  average  radius  of  roll,  which  can  be  expressed  approximately  by  the 
formulas 

S  =  L (1 . 7d  +  Cb*B)  , 

1  |  S 

rs  =  ti'  )  (0.887  +  0.145Cb)e-  2  OG 

The  symbols  L  ,  d  ,  B  and  Ca  represent  length,  draft,  beam  and  block 
coefficient  of  the  ship,  respectively.  The  vertical  distance  from  the  origin 
to  the  center  of  gravity,  OG  ,  is  measured  downward  as  positive. 

Since  in  Kato’s  formula  the  friction  damping  coefficient  does  not  include 
the  amplitude  in  the  case  of  laminar  flow,  the  damping  is  apparently 

linear.  This  is  because  the  damping  is  originally  defined  as  a  nonlinear  form 
and  then  the  nonlinearity  is  cancelled  by  the  dependence  of  the  coefficient 
on  the  roll  Reynolds  number. 

Recently  Myrhaug  and  Sand  [32]  have  carried  out  the  boundary  layer  calcu¬ 
lation  on  a  rolling  circular  cylinder  and  obtained  almost  the  same  result  as 
Kato  by  solving  a  Stokes-type  equation  in  case  of  laminar  flow.  Iheda  et  al . 

[33]  have  also  found  the  same  relationship  between  the  Stoki  "elution  and 

the  Blasius  formula  applied  to  the  unsteady  flow.  Further,  thev  have  confirmed 
the  validity  of  Kato's  formula  in  practical  use,  through  the  measurements  of 
the  velocity  profile  in  the  boundary  layer  on  two-dimensional  cylinders  of 
shiplike  sections. 

In  the  presence  of  forward  speed,  Schmitke  [31]  has  applied  the  skin- 
friction  law  of  turbulent  flow  directly  to  the  case  of  ship  rolling,  in  a 
manner  similar  to  Kato's  treatment.  However,  we  here  cite  Tamiya's  formula 

[34]  based  on  the  more  rigorous  analysis  of  three-dimensional  boundary  layers 


(4.4) 

(4.5) 
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on  rolling  cylinder.  It  can  bs  expressed  in  the  following  simple  form 


BP0  (1 


,4l  JL) 

coL  1 


(4.6) 


where  the  constant  4.1  has  been  determined  through  experiments  on  elongated 
spheroids  in  roll  motion.  The  coefficient  represents  the  friction  damp¬ 

ing  at  zero  fox-ward  speed,  which  can  be  predicted  by  Kato's  formula. 

The  Tamiya  formula  has  been  confirmed  as  giving  a  good  prediction  by 
Ikeda  et  al.  [35],  who  carried  out  the  somewhat  detailed  calculation  of  the 
three-dimensional  turbulent  boundary  layer  development  on  the  axisymmetric 
body  in  its  roll  motion.  Fig.  4.1  shows  a  comparison  of  the  frictional  damping 
as  given  by  the  prediction  methods  and  by  measurements  on  a  spheroid  model.  He 
can  observe  that  the  Kato-Tamiya  prediction  method  is  in  reasonable  agreement 
with  the  measurements. 
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Figure  4.1.  Frictional  component  of  roll  damping  force 


It  has  been  claimed  that  this  Kato-Tamiya  method,  Eq.  (4.3)  through  Eq. 
(4.6),  can  safely  be  applied  to  the  case  of  an  actual  ship  hull  form,  since 
the  ratio  of  friction  damping  to  total  damping  is  usually  quite  small.  The 
exact  treatment  may  be  impossible  theoretically  and  experimentally. 
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According  to  the  above  formula,  the  increment  of  frictional  damping  due 

to  snip  forward  speed  is  proportional  to  the  speed.  When  the  frequency  tends 

to  zero  at  a  finite  speed,  however,  the  B_  value  becomes  infinity,  because 

■fc 

of  the  dependence  of  the  value  B  on  the  roll  Reynolds  number.  This  point 

r  Q 

may  be  a  problem  to  be  studied. 

It  is  noted  that  the  scale  effect  is  quite  large,  since,  at  the  correspond¬ 
ing  speed  and  the  frequency,  the  non-dimensional  value  of  the  friction  damping 
is  almost  inversely  proportional  to  X0,75  ,  where  X  represents  the  scale 
ratio  of  the  ship  to  the  model.  Therefore  the  value  for  the  actual  ship  becomes 
1/20  v  1/30  times  that  for  the  model,  so  that  we  can  safely  ignore  the  friction 
damping  of  the  actual  ship. 

4.3  Eddy  Damping 

In  the  absence  of  ship  speed,  this  component  is  caused  by  flow  separation 
at  the  bottom  of  the  ship  hull  near  stem  and  stem  or  at  the  bilge  circle  near 
the  midship  portion.  The  pressure  drop  in  the  separation  region  gives  rise  to 
this  damping.  Since  early  times,  this  has  been  treated  in  a  manner  similar  to 
that  for  the  drag  problem  in  steady  flow.  Works  by  Bertin  [25]  ,  Watanabe  and 
Inoue  [11] ,  and  Tanaka  [  4  ]  have  been  carried  out  in  this  manner,  in  which  the 
damping  is  assumed  to  have  nonlinear  form  B2U)i  ,  with  the  coefficient  b2 
constant,  depending  only  on  the  ship  hull  configuration. 

In  recent  times,  however,  it  has  been  found  that  the  drag  coefficient  of 
a  body  in  an  oscillatory  motion  varies  with  the  amplitude  of  the  oscillation. 

The  same  situation  may  occur  in  the  case  of  roll  damping.  Ikeda  et  al.  [36] 
investigated  this  point  experimentally  for  a  number  of  two-dimensional  cylinders 
with  ship-like  sections.  In  the  experiments,  the  eddy  dancing  was  obtained  by 
subtracting  from  the  total  measured  damping  (i)  the  wave  damping,  calculated  or 
derived  from  the  measured  radiation  wave  height,  and  (ii)  the  calculated 
friction  damping. 

They  confirmed  through  the  analysis  that  the  eddy  damping  coefficient  can 
safely  be  considered  to  be  constant  in  case  of  ship  rolling.  Figs.  4.2  and 

4.3  show  a  couple  of  examples  of  the  damping  coefficient  BgQ 
two-dimensional  cylinders  of  shiplike  sections. 


measured  for 
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They  further  proposed  a  new  formula  for  the  eddy  damping  for  ordinary  ship 
hull  forms  in  a  kind  of  extension  of  Watanabe  and  Inoue's  approach.  This  can 
be  written  in  terms  of  the  2-D  cross-sectional  coefficient: 


se0 


4  >  . 

~  Pd  <»4»a 


frmax]  2  f R 

[an-  V 


Ho  ,o 


OG 


°P 


(4.7) 


where  d,  r  ,  R  ,  K  ,  a  and  OG  denote,  respectively,  draft,  the 
max  o 

maximum  distance  from  the  point  G  to  the  hull  surface,  bilge  radius,  half  the 
beam-draft  ratio,  area  coefficient  of  the  section,  and  distance  from  the  point 
0  to  G  (downward  positive) .  Thus  the  function  F  can  be  determined  only 
by  the  hull  shape  and  the  pressure  coefficient  by  the  maximum  relative- 

velocity  rati.o  on  the  hull;  details  will  be  stated  i,.  the  Appendix.  Integrating 
the  sectional  value  over  the  ship  length,  we  can  obtain  the  eddy  damping  value 
for  the  given  ship  form. 

In  the  presence  of  ship  forward  speed,  on  the  other  hand,  the  separated 
eddies  flow  away  downstream,  with  the  result  that  nonlinear  damping  decreases. 
Actually,  linear  lift  damping  prevails,  as  can  be  seen  later,  but  we  briefly 
consider  this  decrease  of  eddy  damping. 

Since  the  theoretical  treatment  is  quite  difficult,  we  can  verify  from 
the  results  of  experiments  that  the  amount  of  the  decrease  depends  on  the 
body  shape  as  well  as  on  the  aspect  ratio  of  the  body.  An  example  of  this  is 
show  in  Fig.  4.4,  where  the  abscissa  is  the  reciprocal  of  the  reduced  frequency 
k  =  —  .  In  the  figure,  eddy  damping  at  forward  speed  has  been  derived 

by  subtracting  the  lift  damping  (predicted  separately)  from  the  total  measured 
damping;  the  wave  effect  has  been  excluded  by  covering  the  water  surface  with 
flat  plates.  For  the  case  of  the  ship  form  in  the  figure,  we  can  transfer  the 
value  of  the  extinction  coefficient  b  ,  which  can  be  obtained  by  free  roll 
test,  to  tlie  eddy  damping  coefficient  Bg  ,  because  only  the  eddy  damping 
among  the  component  dampings  is  nonlinear  in  this  case  and  the  coefficient 
does  not  depend  on  the  amplitude. 

In  Fig.  4.4,  the  values  for  the  case  of  the  flat  plate  with  low  aspect 
ratio  represent  the  results  of  the  measurements  for  a  bilge-keel  model  by 
Yuasa  et  ai .  [48],  In  this  case,  when  U  becomes  large,  the  Bg  value  tends 
to  a  constant,  say  B  ,  which  seems  to  depend  on  the  aspect  ratio  of  the  plate. 
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For  blunt  bodies  like  ship  forms,  however,  the  value  B  can  be  regarded  as 

E 

zero.  Xkeda  et  al .  [37]  have  proposed  an  empirical  formula  for  representing 

the  decrease  of  eddy  damping  for  arbitrary  ship  form: 


B 


E  = 


0.04oi2L2 

beq  x  - - — 

lh  +  0.04a/Lz 


(4.8) 


where  the  value  B„ 
Eo 


for  the  zero-speed  case  can  be  predicted  by  Eq.  (4.7). 


Consequently,  the  eddy  damping  for  a  naked  ship  hull  prevails  only  in  the 
absence  of  advance  speed;  it  decreases  rapidly  when  the  ship  moves  forward,  so 
that  it  can  be  neglected  in  the  high-speed  range  of  Fn  >  0.2.  This  fact  has 
also  been  observed  in  the  experiments  by,  say,  Yamanouchi  [12]. 


4.4  Lift  Damping 


Since  a  lift  force  acts  on  the  ship  hull  moving  forward  with  sway  motion, 
we  car:  imagine  that  a  kind  of  lift  effect  occurs  for  ships  in  roll  motion  as 
well.  A  rigorous  treatment  of  this  is,  however,  still  difficult  for  ordinary 
hull  forms. 
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Yumuro  et  al.  [38]  derived  a  simple  formula  by  applying  the  lateral-force 

formula  used  in  the  ship-maneuvering  research  field  to  the  problem  of  roll 

damping.  Accoring  to  his  treatment,  the  damping  moment  M  due  to  lift  effect 

L 

can  be  expressed  in  the  form 


where 


Jtjjj  =  2it_+k(4,1  B/L  -  0.045) 


< 


0  CM  <  0.92 


0.1 

for  0.92  < 

CM < 0.97 

0.3 

0.97  < 

Cfl<  0.99 

(4.9) 

(4.10a) 


(4.10b) 


In  Eq.  (4.9)  and  (4.10),  represents  the  derivative  of  the  lift  coefficient 

of  the  hull  when  towed  obliquely.  The  lever  is  defined  in  such  a  way  that 

the  quantity  corresponds  to  the  incidence  angle  of  the  lifting  body. 

The  other  lever  SL  denotes  the  distance  from  the  point  O  to  the  center  of  lift 
R 

force.  However,  Ikeda  et  al .  [37]  modified  the  values  for  these  levers  origi¬ 
nally  assumed  by  Yumuro  et  al.,  and  they  proposed  another  expression  covering 
cases  when  the  roll  axis  does  not  pass  through  the  point  O  .  The  fined,  pre^ 
diction  form  can  be  expressed  in  terms  of  equivalent  linear  damping: 


(4.11) 


where 


»  0.3d  ,  L  ■  0.5  d 
0  R 


(4.12) 


To  obtain  the  lift-damping  value  experimentally,  one  must  exclude  wave 
effects  by  covering  the  free  surface  with  flat  plates  or  by  carrying  out  the 
measurement  in  a  low  frequency  range  so  that  wave  damping  can  be  neglected. 

Fig.  4.5  shows  a  result  of  measurements  at  low  frequencies,  in  which 
friction  damping  has  been  subtracted  from  the  data  by  using  the  prediction 
formula.  In  this  figure,  therefore,  the  experimental  data  represent  the  sum 
of  the  lift  and  eddy  components.  Since  in  the  high  speed  range  eddy  damping 
can  be  neglected  and,  further,  lift  damping  is  proportional  to  the  advance 
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speed,  we  can  safely  regard  the  data  at  high  speed  as  representing  lift  damping. 
The  solid  line  in  the  figure  shows  the  predicted  values  which  agree  well  with 
the  measured  values  in  spite  of  the  fact  that  the  former  have  been  derived  on 
the  basis  of  a  very  simple  assumption. 


Figure  4.5.  Sum  of  eddy  and  lift  components  of  roll  damping  force 

It  can  be  concluded  that  lift  damping  is  linear  and  that  its  coefficient 
is  independent  of  ai  and  proportional  to  ship  speed,  so  that  it  has  an  impor¬ 
tant  role  in  the  total  damping  at  high  speed.  Particularly  for  ship  forms  such 
as  a  container  ship  or  a  car  ferry,  in  which  the  roll  natural  frequency  is  quite 
low,  lift  damping  becomes  the  most  important  component. 

It  must  be  noted  finally  that  the  prediction  formula  stated  above  may  not 
cover  all  the  varieties  of  ship  forms,  as  pointed  out  by  Ike da  et  ai.  them¬ 
selves.  It  fails  in  cases  of  small  draft-beam  ratio  and  of  ballast  condition  of 
ships,  since  it  is  based  on  the  assumption  that  the  lifting  effect  of  ship  form 
is  approximately  represented  by  that  of  a  flat  plate  with  the  same  length  and 
draft.  It  is  necessary  to  develop  a  more  rigorous  treatment  of  this  component. 
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4.5  Wave  Damping 

In  the  case  of  zero  Froude  number,  the  wave  damping  can  easily  be  obtained 
by  using  the  strip  method,  which  has  been  well  established  recently.  Otherwise 
we  could  even  solve  numerically  the  exact  wave  problem  for  a  three  dimensional 
ship  hull  form.  In  the  strip  method,  the  wave  damping  for  a  ship  section  is 
calculated  from  the  solution  of  a  two-dimensional  wave  problem,  taking  the  form 

*  pNsUw-OG)2  ,  (4.13) 

where  Ns  and  £w  represent  the  sway  damping  coefficient  and  the  moment  lever 
measured  from  the  point  0  due  to  the  sway  damping  force. 

Although  the  direct  measurement  of  wave  roll  damping  is  impossible,  we 
can  use  the  relationship  of  wave  damping  to  the  radiation  wave-amplitude  ratio 
Ar  =  5^/^d .  f;A  the  radiation  wave  amplitude,  and  we  can  compare  it  with  the 
measurements  by  Takaki  and  Tasai  [15].  The  predicted  values  by  strip  theory 
agree  fairly  well  with  toe  measured,  as  shown  in  Fig.  4.6,  There  fore  it  can  be 
considered  that  there  is  little  problem  in  predicting  the  wave  roll  damping  of 

an  ordinary  ship  form  at  zero  speed  by  the  strip  method  with  practical  accu¬ 
racy. 

.  :  calculated 

O  :  measured  at  I.A.M., 


0.0  0.5  1.0  1.5  £d=iu2d/g 

Figure  4.6.  radiation  wave  amplitude  for  Lewis  form  cylinder 


In  the  presence  of  ship  speed,  on  the  other  hand,  it  is  quite  difficult 
to  treat  the  wave  roll  damping  theoretically.  Several  approaches,  for  example, 
Newman  and  Tuck  [39],  Joosen  [40],  Maruo  [41],  Ogilvie  and  Tuck  142],  etc.,  to 
improve  the  strip  method  by  use  of  slender-body  theory  have  recently  been 
attempted.  Troesch  [431  has  evaluated  the  solutions  of  Ogilvie  and  Tuck's 
theory  for  lateral  ship  motions.  However,  a  definite  improvement  in  the  pre¬ 
diction  of  roll  damping  at  forward  speed  has  not  been  reached,  since  Troesch1 s 
result  shows  that  the  first  higher-order  correction  to  roll  damping  is  zero 
according  to  the  slender  body  theory. 

A  method  of  wave-pattern  analysis  to  obtain  the  radiation  potential  of 
ship  motions  at  forward  speed  has  recently  been  proposed  by  Ohkusu  [44].  It 
is  hoped  that  this  method  will  make  it  possible  to  measure  the  roll  wave  damping 
separately.  Otherwise,  we  can  approximate  wave  damping  by  siibtracting  all 
other  predictable  components  from  total  damping  obtained  in  the  forced-roll 
test.  The  results  show  that  wave  damping  behaves  in  a  somewhat  complicated 
manner,  including  hump/hollow  undulations. 

we  can  cite  here  a  couple  of  approximate  treatments  for  predicting  the 
wave  damping  at  forward  speed.  The  first  is  the  method  in  which  the  flow  field 
due  to  roll  motion  is  expressed  by  oscillating  dipoles  with  horizontal  lateral 
axes;  then  roll  damping  i.s  obtained  approximately  from  the  wave-energy  loss  in 
the  far  field..  Hishida  [  9  ]  first  applied  this  treatment  to  the  sway  motion 
of  an  axi symmetric  ellipsoid  to  obtain  the  characteristics  of  wave  roll  damping 
at  forward  speed. 

Ikeda  et  al,  [37]  calculated  the  energy  loss  in  the  far  field  due  to  a 
pair  of  horizontal  doublets  and  compared  the  results  with  experiments  for  models 
of  combined  flat  plates.  Through  these  elementary  analyses  they  proposed  an 
empirical  formula  for  roll  damping  of  ordinary  ship  forms: 


Bw 


®wo 


=  0.  5  [  {  (A.2+  1)  +  (A  -  l)tanh  (  20x  -  b)  } 


+  (2A  -A„  -  1 ) e xp {  -  150  (t  -  0.25)  2}] 


(4.14) 


In  the  second  approximate  method,  the  rolling  ship  hull  is  regarded  as  a 
lifting  body  and  the  wava-wake  problem  for  this  flow  is  solved.  Hanao’ca 
[10]  first  set  the  equation  system  for  the  flat  plate  with  low  aspect  ratio, 
rolling  about  its  longitudinal  axis.  Katanabe  [45]  has  recently  obtained  the 
numerical  solution  of  the  equation,  with  the  result  that  the  roll  damping  is 
expressed  as  a  sum  of  lift  and  wave  damping  in  the  form 

B,  4-  Bh  =  F  (?d,t ,  £j)  ,  (4.16) 
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figure  4.8.  Effect  of  advance  speed  on  wave  and  lift  components 


where  the  detailed  expression  is  omitted  here.  Fig.  4.8  represents  a  result 
of  W.itanabe's  calculation,  which  shows  a  tendency  similar  to  that  of  the  former 
method  in  that  there  is  a  hump  at  t  =  1/4  and  the  value  increases  as  the 
frequency  increases.  This  trend  also  agrees  completely  with  that  of  experi¬ 
ments. 


However,  it  appears  that  there  are  still  some  difficulties  to  be  conquered 
in  both  methods :  In  Hanaoka  and  Watanabe ' s  method,  the  lift  damping  is  based 
on  that  of  a  flat  plate,  not  a  ship  hull,  and  in  the  method  of  Ike da  et  al. 
there  is  a  limitation  in  application  to  ship  forms/  particularly  to  the  case  of 
small  draft-beam  ratio.  It  is  hoped  to  establish  a  more  rigorous  treatment. 

4.6  Bilge-Keel  Damping 

As  stated  in  the  preceding  section,  bilge-keel  damping  is  defined  as  an 
increment  of  damping  when  bilge  keels  are  installed.  It  therefore  includes 
not  only  the  damping  of  the  bilge  keels  themselves  but  also  all  the  inter¬ 
action  effects  among  the  bilge  keels,  the  hull  and  the  waves.  Tanaka  [  4  ] 
and  Kato  [  5  ]  separately  proposed  empirical  formulas  for  bilge-keel  damping 
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in  v  iich  the  effect  of  a  variety  of  ship  farms  was  partly  taken  into  account 
in  terms  of  a  modification  coefficient.  Cox  and  Lloyd  [46]  also  obtained  a 
formula  for  the  bilge-keel  drag  at  zero  ship  speed  using  the  Martin  [  3  ]  and 
Ridjanovic  [47]  experimental  data.  Sasajima  [  6  ]  attempted  to  formulate 
the  hull  pressure  change  due  to  eddy  shedding  from  the  edge  of  a  bilge  keel. 
Watanabe  and  Inoue  [11]  ,  as  mentioned  previously,  also  dealt  with  this  problem 
as  an  extension  of  Bryan's  treatment.  Hishida  [  9  ]  discussed  wave-damping 
due  to  bilge  keels  in  terms  of  a  pair  of  dipoles  on  the  hull  with  axes  tangen¬ 
tial  to  the  hull. 

It  can  be  concluded  through  these  works  that  bilge-keel  damping  is  not 
merely  a  quadratic  nonlinear  form,  but  that  it  depends  on  the  roll  amplitude 
and  the  frequency  in  a  more  complicated  manner,  and  further  that  the  effect 
of  ship  forward  speed  is  not  so  large  as  we  might  expect. 

The  physical  meanings  of  these  facts  have  been  much  more  clarified  by  the 
recent  works  of  Yuasa  et  al.  [48],  Ikeda  et  al .  [49],  and  Fujino  et  al.  [50], 
which  are  based  on  recent  developments  in  research  on  bluff-body  drag  in 
oscillatory  motion.  These  works  will  be  described  in  the  subsequent  sections. 

4.7  Normal-Force  Dap  .ping  of  Bilge  Keel 

To  begin  with,  let  us  consider  the  case  of  zero  ship  speed.  Much  work  on 
the  drag  force  on  a  bluff  body  in  oscillatory  motion  has  recently  been  carried 
out,  mainly  in  the  ocean  engineering  field  [51],  [52],  [53],  [54].  It  has 

also  been  attempted  to  apply  the  results  of  these  works  to  the  problem  of 
bilge-keel  drag. 

Let  the  coefficient  Cp  of  a  body  be  defined  in  the  form 

F  *  CD  .  p  Ajvjv  ,  (4.17) 

where  F  represents  the  drag  force,  ,  A  the  area  of  the  body  projected  onto 
the  crossplane  normal  to  direction  of  motion,  and  v  the  velocity  of  motion. 
Although  CD  is  assumed  to  be  constant  during  the  specified  motion,  its  value 
is  known  to  vary  with  the  period  parameter  or  Keulegan-Carpeter  number,  VT/D 
(V  is  the  maximum  speed  |vmax|  ,  t  is  the  period  and  D  the  maximum  pro¬ 
jected  breadth) ,  contrary  to  the  case  of  steady  flow.  In  the  case  of  periodic 
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oscillatory  motion,  especially  of  bilge-keel  motion,  the  parameter  becomes 
7rr<{>A/bBK  when  we  substitute  r4> auj  ,  kl L  and  2b0K  for  V  ,  T  and  D 
respectively,  where  r  represents  the  mean  distance  from  G  to  the  bilge 
keel  and  b  the  breadth  of  the  bilge  keel-  this  means  that  the  parameter 
can  be  regarded  as  a  sort  of  amplitude  ratio  since  it  depends  no  longer  on 
the  period  of  the  oscillation. 

The  drag  of  the  bilge  keel  can  be  expressed  by  the  following  formula, 
which  was  obtained  by  Ikeda  et  al.  [35]  (including  the  case  of  an  oscillating 
flat  plate),  as  shown  in  Fig.  4.9: 

b 

BK 

Cn  =  22.5  — +2.40  .  (4.18) 

D  tr<{iA 


Figure  4.9.  Drag  coefficient  of  bilge  keel 


Its  validity  has  also  been  confirmed  by  Takaki's  experiment  [54]  except  for 
large  amplitude,  and  a  dependence  on  the  amplitude  ratio  has  appeared  also 
in  the  Cox  and  Lloyd  formula. 

However,  in  order  to  obtain  the  normal-force  damping  of  the  bilge  keel 
installed  on  a  ship  hull  with  comparatively  small  bilge  radius,  it  is  neces¬ 
sary  to  consider  some  modifications  to  Eqs.  (4.17)  and  (4.18).  Ikeda  et  al. 
[34]  assumed  that  the  area  A  in  Eq.  (4.17)  can  be  replaced  by  b 

BK 


per 
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unit  length  and  the  velocity  should  be  multiplied  by  an  empirical  coefficient 
f  of  velocity  increment  at  the  bilge  circle.  Then  the  damping  takes  the  form 


BBKNq 


8 

3tt 


pr2b2  (o 
p  BK 


Iff--- 


W ' 


(4.19) 


f  =  1+0.3  exp  |  -160(1  -  a)! 


(4.20) 


where  a  denotes  the  area  coefficient  of  the  ship  hull  section. 

Fujino  et  al .  [50]  have  recently  assumed  in  a  somewhat  different  way  that 

the  quantity  D  in  the  definition  of  the  period  parameter  should  be  replaced 
RbBK 

by  b  =  b  +  - — -  ,  where  R  is  the  bilge-circle  radius;  they  include 

X  BK  R  +  ^BK 

the  effect  of  the  mirror  image  with  respect  to  the  bilge  circle,  ^instead  of 
taking  D  =  2bBK  '  as  I'cec^a  et  Moreover,  Fujino  et  al.  defined  f  as 

the  average  velocity  increment  at  the  position  of  the  bilge  keel  and  they 
calculated  its  value  for  several  ship  sections  using  the  finite  element  method. 
The  expression  for  damping,  however,  takes  a  form  similar  to  Eq.  (4.19). 

Considering  the  contribution  of  the  amplitude  <j>A  in  Eq.  (4.19)  ,  the 
first  term  in  the  bracket  is  independent  of  the  amplitude ,  so  that  it  becomes 
apparently  a  linear  term.  This  is  because  the  effect  of  period  parameter  on 
the  drag  coefficient  occurs  in  the  damping  terms,  or,  in  other  words,  the 
B2  coefficient  varies  with  amplitude  so  that  a  part  of  B2  is  seemingly  trans¬ 
ferred  into  the  coefficient.  Therefore  Eq.  (2.17),  expressing  the  term- 
by-term  comparison  between  the  damping  and  extinction  coefficients,  does  not 
hold  in  this  case  in  the  sense  of  their  original  definition.  For  the  case  of 
a  naked  hull,  as  seen  previously,  the  effect  of  period  parameter  is  fortunately 
snail,  so  that  we  can  safely  take  B2  as  a  constant . 

The  values  predicted  by  Eq.  (4.19)  for  the  case  of  a  circular  cylinder 
with  bilge  keels  agrees  fairly  well  with  the  measured  values,  as  shown  in  Fig. 
4.10. 


In  the  presence  of  ship  forward  speed,  on  the  other  hand,  it  is  known 
from  the  work  of  Yuasa  et  al .  [48]  that  the  CD  value  of  bilge  keels 
decreases  slightly,  as  can  be  seen  in  Fig.  4.4.  Instead,  linear  damping 


appears,  due  to  the  lifting  effect  of  bilge  keels.  Yuasa  et  al .  applied  the 
approximation  of  low-aspect-ratio  wings  to  the  case  of  bilge  keels  with  the 
result  that  the  normal-force  damping  can  be  expressed  in  the  form 


Figure  4.10.  Component  due  to  normal  force  on  bilge  keel 

A  similar  form  was  also  adopted  by  Schmitke  [31] .  Although  the  decrease  of 
nonlinear  damping  is  not  represented  in  Eg.  (4.21),  the  contribution  of  the 
linear  lifting  term  seems  to  agree  with  the  experimental  results,  as  shown  in 
Fig.  4.11,  in  which  the  ordinate  Cpg  denotes  an  equivalent  nonlinear  drag 
coefficient  of  bilge  keel. 

Recently  Fujino  et  al .  [50]  have  shown  an  impressive  way  of  numerically 
treating  the  drag  of  bilge  keels  at  forward  speed.  They  have  assumed  that  the 
flow  around  a  bilge  keel  corresponds  to  the  steady  flow  around  a  low-aspect- 
ratio  wing  with  breadth  b^.  as  defined  before  and.  with  length  equal  to  that 
of  the  bilge  keel.  Taking  the  effect  of  the  velocity  increment  at  the  bilge 
circle  into  the  expressions  of  the  Incident  angle  and  the  incoming  velocity 
of  the  wing,  they  have  obtained  the  numerical  solution  of  the  bilge  keel  drag 
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on  the  basis  of  Bollay's  nonlinear  wing  theory  [56].  As  we  can  see  in  Fig. 
4.12,  the  results  seem  to  be  quite  promising.  In  the  figure,  the  estimated 
value  at  zero  ship  speed  comes  from  their  method,  mentioned  previously. 

In  conclusion,  it  can  be  noted  here  that  the  slight  increase  of  bilge 
keel  drag  in  the  high-speed  range  is  due  to  the  linear  or  nonlinear  lifting 
effect;  the  problem  of  the  slight  decrease  of  drag  at  low  speed  is  still  left 
for  future  study.  As  a  practical  treatment  at  the  present  time,  one  might 
assume  that  the  change  of  bilge-keel  drag  with  ship  speed  is  negligible,  since 
these  two  effects  cancel  each  other  for  ordinary  dimensions  of  bilge  keels, 
as  seen  in  Fig.  4.11. 

4.8  Hull-Pressure  Damping  Due  to  Bilge  Keels 

Although  this  component  may  be  a  part  of  the  eddy  damping  due  to  the 
interaction  between  bilge  keels  and  ship  hull,  it  is  convenient  to  consider 
it  as  am  independent  component  for  obtaining  a  prediction  formula.  According 
to  the  works  of  Ikeda  et  al.,  the  difference  in  the  hull  pressure  with  and 
without  bilge  keels  is  regarded  as  an  effect  of  the  bilge  keels.  The  hull 
pressure  difference  p  can  be  defined  as  follows,  by  anolugy  with  Cq 
for  the  case  of  normal-force  damping  of  bilge  keels: 

p  =  Cp  j  p|v^|v^  .  (4.22) 

Here,  represents  the  instantaneous  relative  velocity  at  bilge  circle, 

which  can  be  assumed  to  be  =  fr$  ,  with  the  coefficient  f  standing  for 

the  velocity  increment  ratio,  as  stated  in  the  previous  section.  It  is 

noted  that  the  expression  for  pressure,  Eq.  (4.22),  does  not  include  the 
3$ 

term  p—  ($  is  the  velocity  potential),  since  this  term  seems  at  present 

dt 

to  make  a  small  contribution  to  the  damping  in  the  absence  of  a  free  surface; 
further  study  on  this  point  may  be  needed. 

The  distribution  of  the  pressure-dif ference  coefficient  Cp  on  the 
ship  hull,  Eq.  (4.22),  takes  the  shape  shown  in  Fig.  4.13,  according  to  the 
experiments  by  Ikeda  et  al.  [49]  and  Goda  and  Miyamoto  [57].  The  positive 
pressure  c*  in  front  of  the  bilge  keel  is  not  affected  by  the  displacement 
of  the  bilge-keel  motion,  while  the  negative  value  CU  at  the  rear  of  the 
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bilge  keel  and  its  effective  extent  S0  around  the  hull  depend  strongly  on 
the  value  of  period  parameter,  as  does  the  drag  coefficient  of  the  bilge  keels. 

Ikeda  et  al .  made  assumptions  on  the  shape  of  the  Cp  distribution  and 
on  the  values  of  C*  and  C~  to  establish  a  formula  for  the  damping  due  to 
this  pressure  difference,  as  follows: 


3BKH„ 


3tt 


pr2d2w4> .  f2  I 


1  '  d  2~J 


^  1  0 


d  s 


(4.23) 


B=0.2Sm,  1=0  .  lm,  a -eas  -oef  £  .  =0. 935 


The.  Integration  of  I  must  be  done  around  the  whole  girth,  with  the ‘ integrand 
Cp  multiplied  by  the  moment  lever  about  the  rotation  axis.  The  detailed 

expression  will  be  stated  in  the  Appendix.  Ihe  value  of  I  can  be  determined 
as  a  function  of  the  hull  shape  and  the  period  parameter  of  the  motion.  There¬ 
fore  it  includes  not  only  the  apparent  nonlinear  term  of  /  similar  to  that 

of  the  previous  section,  but  also  the  higher-order  terms  of  the  third  and 
fourth  powers  of  ,  although  their  contribution  is  small. 
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Since  it  is  difficult  to  measure  this  component  separately,  the  prediction 
for  the  total  bilge-keel  damping  without  the  effects  of  surface  waves,  that  is, 
BBkjj  plus  ^bkh  '  can  k*2  comPared  with  experiments,  as  shown  in  Fig.  4.14. 

The  agreement  between  them  seems  to  be  fairly  good. 


Figure  4.14.  Effect  of  bilge  keel  on  roll  damping  coefficient 

at  zero  Froude  number 

The  research  for  this  component  in  the  presence  of  ship  forward  speed 
does  not  seem  to  have  been  done.  However,  by  analogy  with  the  case  of  the 
drag  coefficient  one  may  suppose  the  pressure  difference  between  front  and 
rear  surfaces  of  a  bilge  keel  decreases  slightly  and  that  the  width  SQ  of 
its  distribution  on  the  hull  also  decreases  slightly  with  forward  speed. 
Consequently  hull -pressure  damping  decreases  by  a  corresponding  amount.  This 
amount,  however,  will  compensate  the  increment  of  normal -force  damping  due  to 
toe  lifting  effect,  since  it  is  known  that  total  damping  of  bilge  keels  does 
not  vary  much.  We  have  just  to  remember  here  that  the  nature  of  bilge-keel 
damping  changes  frcoi  being  nonlinear  to  being  more  nearly  linear  as  ship  speed 
increases. 
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4.9  Wave  Damping  of  Bilge  Keel 

As  an  example  of  the  study  of  this  component  in  early  times,  we  cite 
Hishida's  theoretical  treatment  [9].  bet  the  radiation  waves  due  to  the 
roll  motion  of  a  bare  two-dimensional  cylinder  be  described  by  C^cos (k0y-uit) 
and  the  change  when  bilge  keels  are  .installed  by  C  ..cos  (k0y-«jt+e )  ,  where 
k^  represents  the  wave  number,  w  /g  ,  y  is  the  horizontal  axis,  t  the 
time  and  e  a  phase  lag.  Then  the  amplitude  of  the  resultant  wave  can 

be  expressed  in  the  following  form  : 


C2  =  c2 

A 


+  SBK  +  2SH  CEKC°Se 


(4.24) 


According  to  Hishida's  consclusion,  although  the  value  of  £  may  be  as 

BK 

large  as  that  of  in  some  cases,  the  magnitude  of  5^  and  thus  of  the 

damping  may  not  always  increase  with  the  addition  of  bilge  keels,  because  of 
the  difference  of  phase  of  the  wave  systems. 

Takaki  [23]  has  recently  carried  out  the  calculation  of  roll  damping  of 
ship-like  sections  with  bilge  keels  vising  the  close-tit  method,  from  the 
results,  as  shown  in  Fig.  4. IS,  it  can  be  said  that  linear  theory  agrees  with 
experiments  only  in  case  of  small  roll  amplitude.  There  is  still  a  consider¬ 
able  discrepancy  for  larger  amplitude,  probably  due  to  the  non-linear  effect 
of  wave  damping. 

However,  it  can  be  noted  that  for  bilge  keels  with  ordinary  breadth  of 
^BK  “  b/6Q  to  B/80  ,  we  can  safely  neglect  the  wave  effect  of  bilge  keels , 
since  the  contribution  of  this  component  is  usually  quite  small  compared  to 
the  viscous  damping  caused  by  bilge  keels. 

In  the  presence  of  ship  speed,  no  study  of  this  component  has  been  made. 
For  ship  forms  with  relatively  large-sized  bilge  keels,  like  a  warship,  it 
may  be  necessary  to  take  the  wave  effect  into  consideration  in  cases  both 
with  and  without  advance  speed- 

At  the  present  time,  however,  it  can  be  concluded  that  the  total  damping 
increment  due  to  bilge  keels  can  be  predicted  as  the  sum  of  two  bilge-keel 
dampings ,  ^bkn  ^hd  * 


Figure  4.15.  Radiation  wave  amplitude  for  cylinder 
with  bilge  keels 

4.10  Prediction  of  Total  Damping 

The  component  dampings  that  were  introduced  in  the  preceding  sections 
can  be  summed  up  to  give  the  total,  damping  of  ship  rolling.  Either  Eq.  (4.1) 
or  Eq.  (4.2)  can  be  used: 

Be  =  Ep  +  BE  +  +  Bw  +  bBKN  +  bBKH  +  BBKW  '  (4-l) 

Be  =  Bf  ...  Be  +  Bl  +  B„  +  Bbk  .  (4.2) 

In  addition,  we  can  add  the  contribution  of  appendages  to  these  expressions, 
as  Schmitke  [31]  has  done.  Although  not  every  individual  component  has  been 
well-established  for  all  ship  forms,  it  is  worthwhile  to  note  that  the  con¬ 
cept  of  dividing  the  damping  into  these  components  should  enable  us  to  predict 
ship  roll  damping  more  accurately  than  by  the  simple  methods.  At  present,  we 
can  indeed  choose  valid  prediction  formulas  for  the  component  dampings;  these 
will  be  identified  in  the  next  section.  Here  we  discuss  the  general  properties 


of  the  prediction  procedure  and  the  component  dampings. 
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In  the  absence  of  ship  speed,  since  only  the  lift  term  B  becomes  zero, 

Li 

the  otner  terms  are  predicted  by  integrating  sectional  values  along  the  length 
of  the  ship,  as  in  the  strip  method  (however,  in  Kato's  formula,  the  friction 
damping  for  the  whole  ship  form  is  still  available).  Therefore  we  can  also 
obtain  the  longitudinal  distribution  of  the  roll  damping.  We  can  compare  it 
with  the  experiments  on  a  divided  model,  as  performed  in  the  Ship  Research 
Institute  of  Japan  [2  3]. 

At  forward  speed,  however,  the  modification  to  the  individual  components 
is  made  for  the  whole  ship  form,  not  for  each  ship  section.  In  this  case  we 
cannot  obtain  the  damping  distribution,  which  presents  an  important  problem 
for  future  study. 

Comparisons  of  these  component  dampings  are  shown  schematically  in  Figure 
4.16  with  the  abscissa  Fn  ,  in  Fig.  4.17  with  frequency  w  ,  and  in  Fig.  4.18 
with  the  amplitude  <p  .  Although  bilge-keel  dampings  are  shown  as  nearly 
constant  in  Fig.  4.16,  the  linear  part  increases  as  Fn  becomes  large.  This 
means  that  in  E'ig.  4.18  the  values  of  bilge  keel  dampings  in  the  limiting  case 
=  0  increase  but  the  slopes  of  the  corresponding  curves  decrease  as  Fn 
increases. 

It  can  also  be  seen  from  Fig.  4.18  that,  the  linear  dampings  consist 

°f  BL  '  BF  '  BW  and  BBKW  '  aS  as  Parts  of  ^KN  and  BBKH  *  The 

nonlinear  part  with  the  constant  B ^  coefficient  includes  be  and  parts  of 

BBK  *  ^ erm3  bp  and  also  depend  on  the  frequency  as  shown  in 

Fig.  4.17,  although  they  appear  to  be  linear  with  respect  to  amplitude.  This 
can  be  realized  from  the  facts  that,  they  were  originally  defined  as  being 
nonlinear  but  then  the  coefficients  varied  with  Reynolds  number  in  the  case 
of  and  with  the  period  parameter  in  the  case  of  BgK  . 

The  scale  effect  on  roll  damping  is  an  important  problem,  no  detailed 
experiments  on  actual  ships  have  been  carried  out  to  determine  the  component 
dampings,  so  that  we  do  not  yet  have  an  exact  picture.  However  it  can  be 
considered  as  in  the  case  of  an  oscillatory  flat  plate  that  dampings  like 
Be  and  Bbk  that  are  generated  by  separating  bubbles  may  not  be  affected 
by  the  scale  difference  between  model  and  ship.  The  earnest  linear  dampings 
like  Bl  and  Bw  are  of  course  independent  of  the  scale  effect.  Therefore 
only  the  frictional  damping  suffers  from  the  scale  effect.  For  instance, 
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Figure  4.18.  Nonlinear  effect  of  roll  damping  components 

when  the  ship  roll  damping  is  to  be  predicted  by  scaling  up  from,  the  model 
value,  the  B0  value  for  the  ship  can  be  obtained  only  by  replacing  the 
Bf  value  by  that  of  the  actual  ship  under  the  same  conditions  of  Fn  ,  w 
and  (fiA  as  for  the  model,  although  Bp  may  be  neglected  in  the  ship  case. 

Concerning  the  independence  among  component  dampings,  there  still  remain 
problems.  All  the  interaction  effects  have  been  included  in  proper  components 
by  their  definitions  in  Sec.  4.1.  However,  in  the  present  prediction  formulas, 
the  interactions  are  usually  ignored,  except  for  the  hull-pressure  change  due 
to  bilge  keels,  which  can  be  regarded  as  the  largest  ’interaction  term.  'Xhe 
other  interaction  effects  may  exist,  and  3Q  they  may  affect  the  independence 
of  components.  However,  we  can  expect  their  magnitudes  not  to  be  large,  since 
almost  all  the  component  dampings  have  been  determined  and  formulated  through 
corresponding  experiments  that  already  implied  the  interaction  effects. 
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4.11  Comparison  with  Experiment 

Although  it  is  still  necessary  to  improve  the  prediction  method  for 
component  dampings,  we  can  at  present  adopt  an  adequate  set  of  prediction 
formulas  for  component  dampings  for  ordinary  ship  hull  forms.  The  following 
is  an  example  of  the  second  (best)  method  which  can  be  chosen  in  a  practical 
point  of  view. 

At  first  we  can  safely  use  Kato  and  Tamiya's  formula  for  the  frictional 
damping  Bp  .  Ikeda  and  ethers'  method  can  be  taken  for  .  The  lift 

damping  BL  is  given  by  the  formula  of  Yumuro,  modified  by  Ikeda  and  others. 

The  wave  damping  can  be  obtained  by  the  strip  method  in  case  of  zero  Froude 
number  and  by  Ikeda  and  others'  formula  at  advance  speed.  The  latter  is  pre¬ 
ferred  over  Hanaoka  and  Watanabe's  method,  which  does  not  yet  yield  an 
adequate  result  in  magnitude,  in  spite  of  the  fact  that  it  is  based  on  a  more 
rigorous  assumption. 

For  the  bilge  keel  damping,  Fujino1 s  treatment  seems  to  be  the  mo3t 
reasonable  one  at  present  time.  However,  since,  it  does  not  Include  hull- 
pressure  damping  due  to  bilge  keels,  B„._.  ,  an  alternative  way  is  to  adopt 

HKxi 

Ikeda  and  others'  method  for  predicting  EBKM  and  BBXH  at  zero  forward  speed 
and  to  assume  that  the  bilge  keel  damping  does  not  vary  with  forward  speed. 
Furthermore  we  can  neglect  both  the  wave  damping  of  bilge  keels  Bbkjj  and  the 
effect  of  the  other  appendages  on  the  hull. 

Using  the  set  of  formulas  thus  chosen,  we  can  calculate  the  roll  damping 
of  a  ship  once  the  hull  form  and  the  roll  conditions  like  Fr,  ,  oo  ,  and  <f‘A 
are  given.  A  computer  program  for  this  is  presented  in  the  Appendix. 

Examples  of  results  are  shown  in  Figs.  4.19  and  4.20  with  the  abscissa 
Fn  and  in  Figs.  4.21  and  4.22  with  the  frequency  as  the  abscissa.  We  can 
recognize  from  FigB.  4.19,  4.20  and  4.21  that  the  magnitude  and  even  toe 
tendency  of  the  roll  damping  are  predicted  by  this  method  to  sufficient 
practical  accuracy.  At  the  same  time,  however.  Fig.  4.22  shows  that  there 
still  appears  some  discrepancy  in  the  trend,  probably  of  the  wave  damping 
at  forward  speed.  Moreover  it  has  been  reported  by  Ikeda  et  aZ.  themselves 
[37]  that  this  method  fails  in  cases  of  small  draft-beam  ratio  or  in  ballast 
condition,  because  of  the  difficulty  of  predicting  the  lift  damping. 


5. 


TREATMENT  OF  NONLINEAR  ROLL  DAMPING  IN  PREDICTION  OF  ROLL  MOTION 


Since  "oil  damping  of  a  ship  is  quite  nonlinear,  as  we  have  seen  in  the 
preceding  chapters,  special  treatment  is  necessary  in  the  calculation  of  ship 
roll  response  in  regular  or  irregular  seas.  Although  there  will  appear  other 
similar  nonlinearities  in  the  restoring  moment  and  the  wave  excitation,  we 
concern  ourselves  hereafter  only  with  the  nonlinearity  of  roll  damping  .n  the 
analysis  of  roll  motion. 


To  begin  with,  the  roll  response  in  regular  seas  can  usually  be  obtained 
from  solution  of  the  simultaneous  equations  of  ship  sway,  roll  and  yaw  motions, 
for  example,  in  the  following  linear  form, 
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where  y  ,  cp  and  tj i  denote  the  displacement  in  sway,  roll  and  yaw,  respec¬ 
tively,  all  expressed  in  complex  form.  In  Eq.  (5.1)  the  coefficients  of  wave 
exciting  forces,  F2  ,  ,  Mg ,  are  expressed  by  an  appropriate  method  as 

functions  of  ship  form  and  speed,  the  wave  amplitude  hA  ,  encounter  angle, 
and  encounter  frequency  we  .  The  inertial  coefficients  in  the  radiation 

forces,  A  • .  (i,j= 2,4,6)  ,  can  also  be  calculated  by  a  similar  method  as 

V 

functions  of  we  .  Then  we  can  easily  solve  the  equations  to  obtain  the 
responses,  provided  that  the  equations  are  linear. 

However  the  roll  damping  coef ficient,  B44  ,  is  usually  expressed  in  a 
form  implying  dependence  on  the  roll  amplitude, 

B44  “  Be(U,  ,  *A).  (5.2) 

In  this  case  the  equation  is  nonlinear'  and  must  be  solved  by  some  kind  of 
iteration  technique.  Then  the  responses  thus  obtained  depend  on  wave  amplitude 
^A  ^  ^  nonlinear  manner.  This  means  that  the  nondi mensional  roll— response 
amplitude  4>AAhA  (k:  the  wave  number)  should  be  stated  for  the  h#  valiae 
concerned. 
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Takaki  and  Tasai  [5]  adopted  the  quadratic  or  cubic  form, 

Be  “  +  ^°^A32  +  T2*??!  '  (2‘6) 

to  express  the  roll  damping  coefficient  Be  of  Eq.  (5.2)  and  they  determined 
the  values  of  and/or  B ^  by  the  experiment.  If  we  use  the  predic¬ 

tion  formulas  stated  above  to  obtain  the  constant  values  of  B1  ,  etc.,  in 
Eq.  (2.6),  however,  a  regression  analysis  with  wide  variations  of  frequency 
and  amplitude  is  needed,  since  the  present  prediction  formulas  are  not  always 
expressed  in  the  simple  form  of  Eq.  (2.6)  but  are  given  as  a  combination  of 
formulas  for  individual  components  in  somewhat  complicated  manner.  In  order 
to  obtain  only  the  roll  response  in  regular  seas,  such  a  transfer  of  damping 
expression  should  not  be  necessary.  We  need  only  the  value  of  the  equivalent 
linear  damping  at  the  specified  values  of  ship  speed,  frequency  and  amplitude, 
as  in  Eq.  (5.2) . 

We  now  proceed  to  the  case  of  irregular  seas.  To  predict  the  ship  response 
in  short-term  irregular  seas,  the  energy-spectrum  analysis  based  on  the  princi¬ 
ple  of  linear  superposition  has  been  applied.  let  the  energy  spectrum  of  the 
incoming  long-crested  waves  be  denoted  by  Sw(cj).  Then  we  obtain  the  energy 
spectrum  of  the  ship  response  5  (a>)  in  the  form 

S(w)  =  |a(u>)  |2SW(W)  .  (5.3) 

The  frequency  cu  of  the  waves  can  be  transformed  into  the  encounter  frequecy 
<oe  of  the  ship.  The  response  amplitude  operator  A(w)  represents  the  ship 
response  per  unit  magnitude  of  input  wave.  In  particular,  it  becomes  <t>A/khA 
(ks  the  wave  number)  in  the  case  of  roll  response. 

The  equation  system  is  nonlinear  when  a(uj)  depends  on  the  amplitude  of 
the  wave  and  thus  on  the  spectrum  S^Cto)  .  Several  attests  have  been  made  to 
treat  this  nonlinearity.  The  most  simple  way,  by  Fukuda  et  al.  [29],  is  to  take 
the  A(io)  value  as  that  in  the  case  of  the  wave  height  (2h^)  equal  to  the 
significant  wave  height,  Hj/3  .  This  corresponds  to  assuming  that  the  roll 
damping  Be  is  constant,  equal  to  the  actual  value  associated  with  a  roll 
amplitude  <j>A  corresponding  to  the  prescribed  wave  height.  The  co-dependence 
of  the  roll  damping,  in  this  case,  is  taken  into  account  exactly. 
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Yamanouchi  [58]  applied  the  regular  perturbation  technique  to  the  non¬ 
linear  terms  in  the  equation  of  motion.  As  another  approach,  the  concept  of 
the  equivalent  linear  damping,  as  stated  in  Chap.  2,  has  been  applied  to  this 
case : 

Be  =  Bi  +V|B2 °4>  ■  <2.13) 


This  has  been  derived  from  the  assumption  that  the  difference  between  linear 
and  nonlinear  damping  can,  in  the  sense  of  least-squares  method,  be  minimized 
during  the  roll  motion,  which  is  assumed  to  be  a  Gaussian  process.  This 
treatment  also  requires  an  iteration  method  of  solution,  in  which  the  variance 
of  the  roll  angular  velocity  is  obtained  in  the  form 
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and  should  be  reevaluated  in  each  iteration  stage.  According  to  Takagi  and 
Yamamoto's  work  [23],  this  method  gives  reasonable  results  in  case  of 
irregular  seas. 


Kaplan  [26]  proposed  that  the  variance  be  replaced  by  the  quantity 


°<J  ~  •  (5.5) 

His  justification  is  that  roll  motion  occurs  largely  in  a  narrow  band  near  the 
natural  freguency  . 

It  should  be  noted  here  again  that  the  coefficients  and  are 

assumed  to  be  constant.  To  apply  this  method,  therefore,  it  is  necessary  to 
rearrange  formulas  to  obtain  a  regression  analysis  on  the  equivalent  linear 
damping  in  a  manner  similar  to  that  stated  before.  For  the  case  of  a  bare 
hull,  however,  this  might  fail,  because  the  roll  damping  is  less  nonlinear  and 
can  hardly  be  expressed  in  a  form  like  Eq.  (2.6)  with  constant  coefficients. 
At  the  same  time,  this  method  cannot  represent  the  w-dependence  of  the  roll 
damping  so  well  as  the  first  simple  approach. 

Dalzell  [59]  ,  Hadd1”  .  [60]  ,  and  Lewison  [61]  have  proposed  expressing  the 
roll  damping  in  a  cubic  torm  without  quadratic  term,  -t- E3$3  ,  for  con¬ 

venience  in  the  statistical  analysis  of  roll  motion  in  irregular  seas. 


6 .  CONCLUSION 


In  th.' s  survey,  recent  trends  and  results  of  research  on  ship  roll 
damping  have  briefly  been  stated.  In  Chapter  2  the  representation  of  roll 
damping  coefficients  was  summarized  and  some  limitations  on  the  definitions 
of  the  coefficients  were  discussed.  The  simple  methods  stated  in  Chapter  3 
could  give  a  rough  estimation  of  the  magnitude  of  roll  damping  at  an  early 
stage  of  ship  design. 

Chapter  4  wa3  cones med  with  the  importance  of  the  concept  of  component 
dampings  and  recent  works  on  them.  Then  a  method  of  predicting  the  shipi  roll 
damping  of  ordinary  hull  forms  was  introduced  and  compared  with  experiments, 
the  resulting  agreement  being  reasonable  from  a  practical  point  of  view. 
Chapter  5  was  devoted  to  comments  on  the  treatment  of  the  nonlinearity  of  roll 
damping  in  the  calculation  of  ship  roll  motion. 

In  conclusion,  it  can  be  said  that  it  is  possible  at  present  to  estimate 
ship  roll  damping  for  ordinary  hull  forms  with  reasonable  accuracy.  However, 
there  still  remain  problems  in  case  of  ballast  condition  of  a  ship  and  for 
unconventional  ship  forms.  From  the  hydrodynamics  point  of  view,  the  linear 
or  nonlinear  lifting  force  on  a  body  in  oscillatory  motion  and  wave  damping 
at  forward  speed  are  still  important  problems. 
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APPENDIX:  An  Example  of  Computer  Program  on  Ship  Roll  Damping 


This  appendix  provides  the  detailed  expressions  that  have  been  omitted  in 
the  previous  chapters.  Then  a  computer  program  of  the  proposed  method  stated 
in  Section  4.11  will  be  shown  as  an  example. 

The  friction  damping  3p  according  to  the  Kato-Tamiya  method  can  be 
written  in  the  form 


-60- 


with 


M  = 


HO 


2(1  +  a  |  +  a^) 


Ho  1  -  OG/d 


a  -  OG/d 
1  -  OG/d 


H  =  1  +  a2  +  9a2  +  2a  (1  -  3a  )cos2t{)  -  6a,  cos  441, 

13  13  ^ 

A  =  -2a3cos  5^  +  a3  (1  -  a3 )  cos  3t(j  +  {  (6  -  3a3  )  a3  +  (a3  -  3a1 )  a3  +  a  2}cos/ , 
B  -  -2a3  sin  5t|i  +  a^  (1  -  a3 )  sin  34»  +  {  (6  +  3a^ )  a2  +  (3a^  +  a2 )  a2  +  a^}  siniji , 
W  =  «D<i  +  at )  sin  iji  -  a3  sin  3iJ»}2  4-  {  (1  -  a^ )  cosi  +  a^cos  3^}^  1/2 


ax  (1  +  a3) 


4a 


3 


-♦1 

(when 

r  (d»  ‘  >  r  (i„ ) )  . 

max'll'-  max  ^2' 

“* 2 

(when 

r  (*1)  <  r  v  (4*2 ) ) » 

max  1  max  * 

(A9) 


f 3  =  1  +  4  exp  {- 1.65  x  10s  (1  -  a)2}  • 


(A1C) 


The  forward-speed  effect  is  stated  in  Eq.  (4.8). 

The  lift  damping  is  given  in  Eqs.  (4.11)  and  (4.12). 

The  wave  damping  at  zero  speed  should  be  predicted  by  the  strip  method; 
the  computer  program  presented  here  does  not  include  this  term.  The  modifi¬ 
cation  of  wave  damping  to  account  for  forward  speed  is  stated  in  Eqs.  (4.14) 
and  (4.15) . 

For  the  bilge -keel  damping  B0K  ,  only  the  terms  BBKN  and  B0kh  are 

used  here;  we  neglect  the  wave  effect  B  Ihe  assumption  is  made  that 

B  is  constant  with  ship  speed.  The  expression  of  B  is  given  in  Eqs. 
BK  BKN  * 

(4.19)  and  (4.20).  B  can  be  expressed; 

BKH 

B3KH0  =  Wpr2dl+0J^A  *  1  *  f2  ' 


(All) 
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where 


X  = 


i_p!,Qds - A  +  B  Cp  ; 


(A12) 


the  coefficient  f  is  given  in  Eq.  (4.20)  ,  arid 

~+  =  1.2 

b 


°P 


c;  -  C^-Cjj=-22.5  - - —  -1.2 


it  r  f 


■A 


(A13) 


Further, 

A  =  (m3  +  m4  )  mg  -  my  , 


B  = 


(1  -  m^2  (2m3  -  m2  ) 


3(H0  -  0.215m!)  *  6(1 -0.215mi )  +  ml  (m3m5  +  m4n6  > 


mv  -  R/d,  m2  =  OG/d,  m3  =  1  -  mi  -  m2  10^  =  H0  -  , 

0.414H0  +  0.0651m^  -  (0.3S2H0  +0.0106)31! 


(A14) 


mc 


(Hq  -0.215m!)  (1  -0.215m!) 

0.414hq  +  0.0651m2  -  (0.382  +0.0106H0)mi 
“  (Ho  -  0.215mi)  (1  -0.215m!) 

S0/d  -  0.25-nmi  (S0  >  0.251TR)  , 


ml  = 


1  0 

(So  <  0.25tiR)  , 

1  m7  +  0.414mj 

(Sq  >  C, 25ttR)  , 

|n>7+  v'z  {1  -  costSQ/R)}  mj  (sQ  <  0.25ttR) 

-circle  radius 

R  and  the  distance  r 

.  are  given  in 

the  form 

iW-0'”11 

1  IT  -  4 

l K' C  f  D/*  /  f 

1 a 

(h0>i,  p/a>i)  , 

'  B/2 

(Hg<l,  R/d>H0)  , 

(A15) 


-62- 


<A16) 


A  computer  program  can  be  written  ir.  the  following  form: 

(1)  INPUT 

1}  Title  (70  words)  In  put  type  70A1 

2)  L,B,D,NABLA,CB,CM  6F10.0 

L  ;  length  b  ;  breadth  D  ;  draft  (m) 

NABLA  ;  displacement  volume  (m3) 

CB  ;  block  coeft.  CM  ;  midship  coeff. 

3)  NUB  F10.0 

NUE  r  kinematic  viscosity  (m3/s) 

Empirical  formula 

v=0 . 0178/ ( 1+0 . 0336t+0 . 00022 It2 )  cm2/s 
t  ;  temperature  (°C) 

4)  M,N  21  5 

M  ;  the  number  of  calculation  points  of  Fn  (interval: 

Fn=0 .02) 

N  ;  tne  number  of  input  data  of  sections 
(=  the  number  of  data  cards  of  5) 

5)  X (I)  ,HJ2(I)  ,SIG(I)  ,BX(I)  ,DX(I)  5F10.0 

X(T.)  ;  the  number  of  S.S  (0.0  at  A.P.&10.0  at  F.P.) 

H0(l)  ;  half  breadth/draft  (=BX(I)/DX(I) 

SIC  (I )  }  area  coefficient  (=cross  section  area/'BX(I)DX(I) 

BX(I)  ;  breadth  (ra) 

DX(I)  i  draft  (m) 

6)  BBK , XBK1 rSBK2  3F10.0 

BBK  •,  breadth  of  bilge  keel  (m) 

XBK1  ;  SS  number  of  aft  end  of  bilge  keel 
XBK2  ;  SS  number  of  fore  end  of  bilge  keel 

7)  OGD , T , THETA  3F10.0 

CGD  ;  =0G/d 

T  ;  roll  period  (sec) 

THETA  ;  roll  amplitude  (rad) 

8)  BWOSM  ^  F10.0 

EvJOSM  ;  Byf  at  Fn-0.0  calculated  by  potential  theory 

(2)  OUTPUT 

1)  prints  of  input  data 

2)  distribution  of  the  coefficient  CR  of  the  eddy  making  component. 

(Ci^Mp^/ (0 , 5pd4  $|$|  ,  Mjyj  ;  roll  moment  due  to  the  eddy  making 
component ) 

3)  distributions  of  ABbk  and  ABn/ABbk. 

ABbk  «  bBK  f°r  unit  length 
ABjj  ;  Bn  for  uni  t  length 

4)  each  component 
rates  of  each  component 
total  toll  damping  coefficient 


for  each  Fn 
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C  *** ••ESTIMATION  OF  ROLL  DAMPING***** 

C  CODED  8V  Y.IKEDEA. 


REAL  L.NUE.NABLA 

DIMENSION  X(25) ,H0( 25 ) . S IG( 25 ) ,BX(25) ,DX(25) ,B44HAT( 100) . 

*  BWHAT( 100) , BEHAT ( 100) , 8LHAT ( 100) 

*  ,FN( 100) ,TITLE{70) ,BFHAT( 10O) 

READ  (5,100)  (TITLEC I ), I* 1 ,70) 

100  FORMAT  ( 70A 1  ) 

C  * • “PRINCIPAL  DIMENSIONS***  L- LENGTH , B-BREADTH . C-DRAFT , NABLA-VOLUM 

C  *«*  CB”BLOCK  "OEFF . ,CM*MIDSHIP  COEFF . , 

READ  (5,101)  L.B.D.NABLA, CB.CM 

101  FORMAT  (GF10.0) 

READ  (5, 103)  NUE 

103  FORMAT  (F10.0) 

C  ***  NUMBERS  OF  FN  AND  DATAS  •** 

REAO  (5,104)  M.N 

104  FORMAT  (215) 

C  *«*  PARAMETERS  OF  EACH  SECTION  *»*  H0-B/20 . SIG-S/BD . BX-BREADTH 
C  DX-DRAFT 

DO  1  I  -  1 , N 

1  READ  (5,105)  X( I ) , H0( I ) , SIG( I ) , 8X( I ) ,QX( I ) 

105  FORMAT  (5F10.0) 

C  *•*  BILGE  KEELS  DATA  «**  8BK-BREA0TH  OF  B.K..XBX1-X  OF  B.K,  END 

C  XBK2 *X  OF  B.K,  END (FOR) 

READ  (5,106)  BBK.XBK1 .XBK2 

106  FORMAT  (3F10.0) 

S  CONTINUE 

C  ***  CONDITION  »««  OGD*0G/O,T*PERIOD,THETA*AMP.  OF  ROLL 

READ  (5,102)  OGO.T, THETA 

102  FORMAT  (3F1O.0) 

IF  (T.LT. 0.000001)  STOP 
C  ***  WAVE  MAKING  COMPONENT  AT  FN-0  *•* 

REAO  (5,103)  BWOSM 

WRITE  (6,200)  (TITIEU). 1-1.70) 

200  FORMAT  ( 1H 1 , /// , 7X ,SH* *  * •* , 2X . 70A 1 . 2X . 5H» *»» • J 
WRITE  (6.300)  L.B.D.NABLA. CB.CM 

300  FORMAT  < 1H  , // , 5X , 8H* ‘DATA** , 2X , 2HL- , FB . 5 , 2X , 2HB* . F8 . 5 , 2X , 2HD- 

*  , F8 . 5 . 2X , 6HNABL A- , F8 . 5 , 2X . 3HCB* , FS . 5 . 2X , 3HCM* , FB . 5 ) 

WRITE  (6,301)  OGO,  T,  THETA  ,  NUE  .BWO.'.M 

301  FORMAT  (  1H  ,  5X  ,  5H0Q./D*  ,  FB  .  3 , 2X  .  2HT*  .  F8 . 3 . 2X  ,  GHTHETA* ,  F8 . 3 , 3HRA0 

*  ,  2X  ,  4HNUE  ■  ,  F  12 . 8 , 2'' .  6HPVQSM*  ,  F  15 . 10) 

WRITE  ( G , 302  )  B3K . XBK 1 . XBK2 

302  FORMAT  ( 1H  . 5X , 4HBBK- . F8 . 5 , 2X , 5HXBK 1  * , FB . 3 , 2X . 5HXBK2* , F8 . 3 ) 

WRITE  (6,304) 

304  FORMAT  ( 1H  . // , 4X , 25H**DA TAS  OF  EACH  SECTION**) 

00  15  I  *  1 , N 

15  WRITE  (6,303)  X( I ) , H0< I ) , SI G( I ) , BX (I ) , DX ( I ) 

303  FQRMAT( 1H  , 4X , 3HSS* , F8 . 3 , 3X , 3HHQ- . F8 . 5, 3X , 6HSIGMA- , F8 . 5 . 3X . 2HB* 

*  , F8 . 5 , 3X , 2H0* , F8 . 5 ) 

C  ••*********'*‘***********************************‘**'-****** 

OMEGA ”6 . 283 1 8/T 
DO  2  I « 1 . M 

2  FN( I ) *0.0+0 , 02 ‘FLOAT ( I -  1  ) 

CALL  FRICT  < L , B . D. CB .NABLA .OGO . NUE , OMEGA , FN , BFHAT , M ) 

.  CALL  WAVE  (L,0, OMEGA, EWOSM.FN.BWHAT.M) 

CALL  LIFT  (L.B.O, CM, NABLA, OGO. FN.BLHAT.M) 

CALL  EDDY  ( X , HO , S IG , BX , DX , 8 , 0 , NABLA , OGO . OMEGA , THETA , FN , 

*  BEHAT. M.N.L) 
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59 

60 
61 
62 

63 

64 
05 
66 

67 

68 

69 

70 
7  1 

72 

73 

74 

75 

76 

77 

78 

79 

80 
31 
82 
83 


IF  ( BBK , LT . 0 . C00000000 1 )  GO  TO  20 

CALL  BK  (X, HO, S IG , BX . DX , B , D , NABLA , OGD , OMEGA, THETA, BBK , 

*  XBK 1 .  XBK2  ,  BBKHAT  ,I1,N,L) 

20  IF  ( BBK. LT. 0.0000000001)  BBKHAT *0.0 

00  3  I  ■  1 ,  M 

3  B44HAT ( I ) *BFHAT ( 1) +BWHAT ( I ) +BLHAT (I ) +BEHAT ( I )+BBKHAT 
WRITE  (6.201) 

201  FORMAT  ( 1H  . /// . 4X . 2HFN . 6X , 5HBFHAT , SX . 5HBWHAT , 5X . 5HBEHAT , 5X . 

*  5HBLHAT , 5X . 6HBBKHAT , 4X , 2H*  * , 3X , 6HBF/B44 , 4X , 6HBW/B44 , 

*  4X . 6HBE/B44 , 4X . 6HBL/B44 . 2X , 7HBBK/B44 , 4X . 2H*»  . 2X . 

»  6HB44HAT, 2X,2H»*  ) 

00  4  1*1, M 

BFF'BFHATt I ) /B44HAT ( I ) 

BWW-BWHAT( I ) /B44HAT ( I ) 

BEE “B  EHAT ( I )/B44HAT( I) 

BLL "BLHAT ( I )/B44HAT( I ) 

BBKK*B8KHAT  /B44HAT( I ) 

4  WR I TE ( 6 , 202 ) FN( I ) , B FHAT ( I ) , BWHAT ( I  ) . BEHAT ( I ) . BLHAT ( I ) .BBKHAT, 

*  BFF , EWW , BEE , 8LL . BBKK , B44HAT ( I ) 

202  FORMAT  (1H  . 2X . F6 . 3 , 2X . F 8 . 5 , 2X . FB . 5 , 2X , F8 . 5 . 2X , F8 . S . 2X , F8 . 5 

*  , JX,2H**.2X,F8.5.2X.F8.5,2X,F8.5,2X,F0.5.2X,F8.S,2X, 

*  2H**. 1X.F8.5, IX , 2H* • ) 

GO  TO  5 

END 


| 


5 

j 


'i 
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83 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
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^«*tt«4*4f*(»«4«t«*tl«**tt******a»«f<**»4*>«*a**«««l»i***«******i*tt4l***tt»t»*« 

SUBROUTINE  LIFT  ( L . B , 0 , CM , NABLA , OGO , FN , 8LHAT , M ) 

C  LIFT  COMPONENT 

C  REF.  V.IKEDA  ET  AL  <JZK.N0.143) 

REAL  L, NABLA, KAPA.KN.LO.LR 
DIMENSION  FN( 100) ,BLHAT( 100) 

IF  (CM. LE. 0.92)  KAPA«00 

IF  ( CM . LE . 0 . 97 . ANO . CM . GT . 0 . 92 )  KAPA-0.1 

IF  (CM. GT, 0  97)  KAPA-0.3 

KN-6 . 28319*0/L+KAPA*(4 . 1*8/L-G.045) 

OG-OGO’D 
L0«0  3*0 
LR*0 • 5*0 
oa  i  I ■  i . m 

1  SLHAT ( I )*L*O*KN*L0*LR*FN( I  ) *0 . 5/ ( NABLA'B*  *2 ) "SORT (0.3*L*B)*( 1 .0- 
*  1 .4*0G/LR+0.7*0G**2/(L0*LR) ) 

RETURN 

ENO 


<PAGE  3> 


/////  FILE; -A  llllt 


>»»  SUBROUTINE  LIFT  «,« 


-67- 


PaGE  4> 


+ - - - --“---1 

I  DATE :09-25-B1 ,14:  1 6  ! 


102 

103 

104 

105 

106 

107 

108 

109 

1 10 
1  1 1 
1 12 
1  13 

114 
1  15 

115 


SUBROUTINE  WAVE  ( L , 0 , OMEGA . BW05M , FN . BWHAT , M ) 

C  WAVE  MAKING  COMPONENT 

C  REF  Y.IKEOA  ET  AL  (UZK . NO . 243 ) 

REAL  L.LOMEGA 

DIMENSION  FN( 1001 ,BWHAT( 100)  . 

GUZA ID «OMEGA*  *2*0/9 , 60665 
A 1 ■ 1 .O+GUZAID*’ ( -  1 . 2 ) *  EXP ( -2 . O'GUZAID) 

A2*0.5+GUZAI0**( -1 ,0)*£y.P(-2 ,0*GU2AID) 

DO  1  I*1,M 

L0MEGA*0MEGA*FN(I)*3QRT(L/9 . 80665) 

1  EWHAT  (  I  )»BWQ5M*0.5*(((A2*-1  ,0)+(A2-1  .0)  *TANH(  20 . 0*  ( LOMEGA-O.  3  )  )  )  * 
*  (2.0*A1-A2-1 .0) *  EXP t" 150.0* (LOMEGA-O. 25)* *2) ) 

RETURN 

END 
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1  17 
118 
1  19 
120 
121 
122 

123 

124 

125 

126 

127 

128 
129 


c. .......... 

SUBROUTINE  FRICT  ( L . B . 0 . CB  .  NA8LA , OGO , NUE , OMEGA , FN , BFHAT , M ) 

C  FRICTIONAL  COMPONENT 

C  REF.  H.KATO  (J2K.NQ.102)  ANO  S.TAMIYA  ET  AL  1UZK.N0.132) 

REAL  NAB LA , NUE , L 
DIMENSION  FN( 100) ,BFHAT( 100) 

SF-L*( 1 . 7+0+C8-S  ) 

RF*( (0.887+0. 145*CB)«< 1 . 7 -D+C8+B) -2 . O’OGO+D )/3 . 14S 
00  1  1-1. M 

1  BFHAT (I)«O.78  +  SF*9F**2  +  S0RT( OMEGA* NUE *B/ t9,6133)/( NABLA+B*1  2 ) 
*  *( 1 .0+4 . 1 *FN( I  ) /OMEGA -SORT ( 9 . 80665/L ) ) 

RETURN 

END 
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130 

131 
1  32 

133 

134 
13b 

136 

137 

138 

139 

140 

14  1 

142 

143 

144 

145 

146 

147 

148 

149 

150 

15  1 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 
■  173 

174 

175 

176 

177 

178 

179 
1B0 
181 
182 
103 
184 

105 

106 
107 


. . . 

SUBROUTINE  EDDY  (X . HO . SIG, BX , DX , B , D , NABLA ,OGD . OMEGA . THETA . FN, 

*  BEHAT.M.N.L) 

C  EDDY  MAKING  COMPONENT 

C  REF.  Y.IKEDA  ET  AL  ( UZK . NO . 1 42 , UZK . NO . 1 43 ) 

REAL  NABLA. LTHETA.L 

DIMENSION  X ( 25 ) ,H0(2S) .516(25) ,BX(25) ,DX(25) ,FN( 100) ,BEHAT( 100) , 

*  RMAX  I  (2 ) . V( 2 ) , CR { 25) , CR 1 ( 25 ) . X 1 <  30 ) 

00  1  0-1, N 

AHO-HO( J)/(1 .O-OGO) 

SIGMA- ( SIG(J)-OGD)/( 1 ,0-QGD) 

E  * ( AHO-  1  . 0 )/ ( AHQ+ 1 .0) 

E2-E--2 

A  *  4 . 0  *  S I GMA - ( 1 .Q-£2)/3 , 1415  +  E2 
0*-A/(A  +  3.0> 

O2*S0RT(O**2-( A-1 ,0)/(A+3.0) ) 

A3  *0+02 
A1-E‘(  1  . 0+A3 ) 

AM-BX(  J)A 1 ,0+A I+A3) '0 .5 
AA1*A1*(  1  ,0+A3)/A3*0.25 
IF(AAI.QT.I.O)  AA  1-1.0 
IF(AA1 .LT.-1 .0)  AA1--I.0 

DO  2  1-1,2 

LTHETA-0.5*ARCOS(AA1 ) 

IF  (I .EO. 1)  LTHETA-0 . 0 

AH- 1 . 0+ A  1 +'2+9 . 0-A3**2+2.0*A1-( 1 .0-3 . 0-A3 ) -COS( 2 .O-LTHETA ) -6 . 0-A3 

*  *COS( 4 , 0*  L THETA ) 

AA-~2.O*A3*C0S<5.0*l-THETA)+A1*(  1  . 0-A3 ) -COS( 3 . O-LTHETA )+(( 6 . 0-3 . 0* 

*  A 1 )*A3**2+( A 1--2-3 .0-A1 ) * A3+A 1  * *2 ) -COS ( LTH6T A ) 
BB*-2.0*A3'SIN(5.0*LTHETA)+A1*( 1 . 0- A3 ) *S I N( 3 . O-LTHETA )+ ( ( 6 . 0+3 . 0- 

*  At )+A3*-2+(3.0*A1+A1*-2) * A3+A 1  * *2 ) • SIN( LTHET A) 
V(I)*2.0+AM*SQRT(AA-*2+8B»-2)/AH 

%  RMAXii I  ) -AM* SORT U ( 1 .0+A1 ) -SIN( LTHFTA ) -A3-S !N( 3 . O-LTHETA ) )*-2 

*  ♦(( t .  Q-A 1 ) -COS ( LTHET  A )+A3 -COS ( 3 .O-LTHETA ) )**2) 

RMAX-RMAX  1  (  1  ) 

VMAX-V(  1 ) 

IF  (RMAX  1  ( 1 ) , LE . RMAX 1(2})  GO  TO  8 
GO  TO  9 

8  RMAX-RMAX  1(2) 

VMAX-V/C  2  ) 

9  CONTINUE 

RMEAN»2.0*DX( d)*( 1 . O-OGO ) *SQRT( AHO-3 IGMA/3 . 1415) 

P  1  -VMAX/RMEAN 
P2-RMAX/RMEAN 
PP3-P1  +  P2 

IF  ( SIGMA. LT. 0.99)  GOTO  20 

GAMMA- ( 1 . 0+4 . 0*EXP( - 1G5000 . 0+ ( 1 . Q- SI GMA ) --2 ) ) -PP3 
GOTO  21 

20  GAMMA-PP3 

21  CP«Q.3*(0.87*EXP(-GAMMA)-4 ,0*EXP(-0. 107-GAMMA )+3 .0) 

F 1*0.5* ( 1 ,0+TANH(2O.0*(SIG( J)-0. 7) ) ) 

F2-0 .  5*  (  1 . O-COS ( 3 , 1415*SIG(J) ))-1 . 5-( 1 . 0-EXP( -5 . 0* ( 1 .O-SIG(d) ) ) ) * 
«  SIN(3. 141S*5IG(J)  )«*2 

R«2.0*DX(G)*SORT(H0(J)»(SIG(J)-1 ,0)/( -0.8584 ) ) 

RD-R/DX  tO) 

IF  (HO(d).LE. I.O.ANO.RD.GE.AHO)  R-O.S-BX(d) 

IF  (HO(d).GT. 1.0.AND.RD.GE, 1 .0)  R-OX(d) 

RD*R/OX(J) 
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188 

1 

CR  1  ( U ) -UMAX* "2/OX (  j)  '*2  *CP*((  1 , 0-F  1  *RD  )  *(  1  ,0-OGD-F  1  *RD)+F2* 

189 

(H0(J)-F1-R0)*"2) 

190 

WRITE  (6.452) 

19  1 

452 

FORMAT  (  1H  ,//,4X,4H****, 3 1HL0NCITUDI NAL  DISTRIBUTION  OF  CR, 

192 

* 

4H**** ) 

193 

00  10  U-1.N 

194 

10 

WRITE  (6.453)  X(0),CR1(J) 

195 

453 

FORMAT  (  IH  . 4X . 3HSS* , F 8 . 5 , 4X , 3HCR- , F8 . 5 ) 

196 

00  3  K-  1 .21 

197 

X 1 ( K ) »0 . 0+0 ,5*FL0AT(K~ 1 ) 

198 

MAX«N 

199 

3 

CALL  H0KAN1  ( X , CR 1 .  25 . MAX . X 1 (K ) , CR ( K ) . DAM . 1 . 0) 

200 

CR(  1  )«1  .  5*( 1 ,0-OGD) 

201 

CR ( 2  1  )-1 .5*< 1 .O-QGO ) 

202 

SAM-0.0 

203 

00  4  K- 1 . 10 

204 

K2-2-K 

205 

K1-K2-1 

206 

K3-K2+ 1 

207 

SAM1-CR(K1)+4.0*CR(K2)i-CR(K3) 

208 

4 

SAM-SAM+SAM1 

209 

CRT -SAM/ 60 . 0 

2  10 

BEHAT ( 1 )-4.0*L*0**4/3 .0/3. 1415*0MEGA*SQRT(B/19.6)/NABLA/B*»2’CRT 

2  1  1 

* 

"THETA 

212 

DO  5  I-2.M 

2  13 

AK-OMEGA/FNl I ) "SQRT( L/9 . 8  ) 

2  14 

5 

BEHAT ( I ) -BEHAT ( 1 ) • (0 . 04 • AK ) * *2/( (0 . 04 *AK ) » *2+ 1 .0) 

215 

RETURN 

216 

END 

\ 
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254 

255 

256 

257 

258 
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SUBROUTINE  BK  ( X , HO , 5 IG , BX , DX , B . D , NABLA , OGD , OMEGA . THETA . BBK , 

*  XBK  1  , XBK2 , BKHAT , M , N , L ) 

OAMPING  DUE  TO  BILGE  KEELS. 

REF.  Y.IKEDA  ET  AL  (K2K.NO.161,  KZK.N0.16S) 

REAL  NABLA.  M 1 , M2 . M3 , M4 . M5 . M6 , M7 , M8 , L 

DIMENSION  X(25),H0(25),StG(25),BX(25),0X(25), BBKHAT (25).XBK(25) 

*  , RATIO! 25 ) 

X8K{ 1 )-XBK1 

XBK( 1 1 )*XBK2 
DO  1  1-2,  10 

1  XBK( I  ) -XBK (1-1 )  +  ( XBK2-XBK 1 ) *0.  1 
MAX-N 


00  2  1*1,11 

CALI.  H0KAN1  (X..HO,  25  .MAX  .  XBK{  I  )  .  HO  1  ,  DAM.  1  .0) 

CALL  H0KAN1  (X.SIG,  25 . MAX , XBK( I ) . SIG 1 , DAM ,1,0) 

CALL  H0KAN1  (X.OX.  25 , MAX , XBK ( I ) . DX 1 , OAM , 1 . O ) 

CALL  HQKAN1  IX, BX,  25 . MAX . XBK ( I) , BX 1 , OAM . 1 , 0) 

R*2 . 0*DX 1 ’SORT (HO  1  * ( S IG 1 -  1 . 0 ) /( -0 . 8585 ) ) 

RD-R/DX1 

IF  (HOI .LE . 1 .0. AND . RD . GE  HOI }  R-0.5*RX1 

IF  (HOI .GT. 1 .Q.AND.RO.GE. 1 .0)  R-0X1 

RD-R/0X1 

F-1 .0*0. 3-EXP ( -160,0* ( t .0-SIG1 ) ) 

RBK-OX I “SORT ((HQ1-0.2929*RD)**2+( 1 . O-OGD-O . 2929-RD ) * *2 ) 
M1*RD 


M2-0G0 

M3* 1 .0-M1-M2 
M4-H01-M1 

M5* ( 0 . 4 1 4 -H0 1 +0 . 065 1 *M 1  *  *  2- (0 . 38 2* H0 1+0 . 0106 )  *M  1 ) / ( ( HO  1 -0 . 2 1 5 

*  *M1)*( 1.0-0. 215*M1)) 

M6* ( 0. 4 14 *H0 1 +0. 065 1*M 1**2- (0.382+0. 0106 -HOI) ’Ml )/( (HO  1 -0 . 2 15 
»  *M1 )*( 1 .0-0.2 I5*M! ) ) 

S0»0.3-‘(3.  14  15’F’RBK‘THETa  )  +  1  ,95‘BBK 
M7“ SO/OX  1-0.25*3. 14 13’MI 
H1-0.25’3. 1415’R 
IF  (S0.LT.R1)  M7-0.0 

M8-M7+0. 4 14*M1 

IF  (S0.LT.R1)  M8-M7+1 ,414*(  1 . 0-COS ( SO/R ) )*M1 

A*  (M3+M4  ) ’M8-M7  ”2 

BB«M4**3/3./(H01-O.215*M1 )*( 1 .0-M1 )  ”2  * ( 2 . 0’M3-M2 )/6 .0/ ( 1 .0-0.219 

*  *M 1 ) +M 1  * ( M3 “M5+M4  *M6 ) 

CPPLAS* 1 .2 

CPMINS--22 .5*BBK/(3.  14 15*RBK*F’THET A ) -  1 . 2 

CO’CPPLAS-CPMINS 

»*♦  «*  2HAT  FOR  UNIT  LENGTH  ”* 

RATI .  )  *RBK*BBK*CD/( RBK*8BK*CD+0 ,5*DX1**2*( -A*CPMIN5+BB ’CPPLAS ) ) 

2  BBKHAT ( I)»8.0*RBK  ”2 ’OMEGA* SORT ( B/ 19 . 6 ) ’THETA’F " *2/ (3.0*3. 1415 

*  *NABLA*B**2)*(RSK*SBK*CD+O.5*DX1*»2*(-A*CPMINS+0B*CPPLAS 

*  )/ 

WRITE  (o.lOO) 

100  FORMAT  ( 1H  , // , 4H**** . 35HL0NGITUD I NAL  DISTRIBUTION  OF  BBKHAT, 

*  4H**»*j 

00  7  1*1,11 

7  WRITE  ( 6 , 101 )  XBK(  I ) , BBKHAT ( I ) , RAT I0( I ) 

101  FORMAT  (  111  .4X.3HSS*  ,  F8 . 5 , 3X , 7HBBKHAT- , F 1 3 . 8 , 3X , 22HN0KMAL  FORCE/ 
•TOTAL  BK* , F 1 3 . 9 ) 

****  BBKHAT  FOR  THREE  DIMENSIONAL  SHIP  FORM  **** 
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275  SAM  *0.0 

276  DO  3  I  «  1 ,5 

277  1 2 *' 2 "  I 

278  SAMI -a8KHAT( 12-1 ) +4 . O+BEKHAT ( 12  )  +  B8KHAT( 12+1 ) 

279  3  SAM*  SAM+SAM 1 

280  BKHAT*  SAM*  (XEH.2-X8K1  )  *0 .  1/3. 0*1.  *0.  1 

281  RETURN 

282  END 
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V 


I 


i 


L 


i 
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203 

C  + 

*  *  *  + 

234 

SUBROUTINE  H0KAN1  (XI, VI. MAX 

285 

c 

LAGRANGE  3  POINTS  INTERPCLAT 

2S6 

DIMENSION  X 1 C  MAX ) , Y 1 ( MAX ) , WX 

287 

Nt  -N- 1 

203 

OQ  10  J«2.N1 

209 

IF  (X . LE .X 1 ( I)  )  GO  TO  1 

290 

10 

CONTINUE 

29  1 

1 

11*1-1 

292 

IF  (X.GT.XMNO)  I1-N-2 

293 

12*11+2 

1.94 

00  20  I  "I  1 ,  i2 

295 

11*1-^1-11 

296 

WX ( I I ) - X 1  (  I  ) 

297 

20 

WY< II)-VI(I) 

298 

IF(M1  NE. 1)  GO  TO  2 

299 

CALL  LAG3(WX.WY,X.T) 

300 

2 

CONTINUE 

301 

1FCM2.NE.1)  RETURN 

302 

Y  X  =0 . 0 

303 

RETURN 

304 

END 

I  OATE : 09 -25 -81. 14: 16 | 
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305 

c*  *  *  ** 

306 

SUBR0UTTNE  LAG3  (WX.WY.X.Y) 

307 

DIMENSION  WX ( 3 ) ,WY(3) 

308 

Y  aQ  .  0 

309 

DO  11  1-1,3 

3  10 

W-  1  .0 

31  1 

Z«1  .0 

312 

DO  12  U-1,3 

313 

IF(O.EO.I)  GO  TO  12 

314 

W«W*(X-WX( J) ) 

315 

Z«ZMWX(  I  )-WX(  J) ) 

316 

12 

CONTINUE 

317 

Y - Y+WY ( I )-W/Z 

318 

1 1 

CONTINUE 

319 

RETURN 

320 

END 

i 
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